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Electron dynamics in the bulk and at the surface of solid materials are well known to play a key 
role in a variety of physical and chemical phenomena. In this article we describe the main aspects 
of the interaction of low-energy electrons with solids, and report extensive calculations of inelastic 
lifetimes of both low-energy electrons in bulk materials and image-potential states at metal surfaces. 
New calculations of inelastic lifetimes in a homogeneous electron gas are presented, by using various 
well-known representations of the electronic response of the medium. Band-structure calculations, 
which have been recently carried out by the authors and collaborators, are reviewed, and future 
work is addressed. 
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I. INTRODUCTION 

Over the years, electron scattering processes in the 
bulk and at the surface of solid materials have been the 



subject of a greatj-|V|ariety of experimental and theoreti- 
cal investigations.ErEl Electron inelastic mean free paths 
(IMFP) and attenuation lengths have been shown to play 
a key role in photpelectron spectroscopy and quantitative 
surface analysis dQ Linewidths of bulk excited electron 
states in metals have also been meatiUEed, with the use 
of inverse photoelectron spectroscopyJjilj More recently, 
with the advent xtfjtime- resolved two-photon photoemis- 
sion (TR-2PPE)E3o and ultrafast laser technology, time 
domain measurements of the lifetimes of photoexcited 
electrons with energies below the vacuum level have been 
performed. In these experimepts^the lifetimes of both 
hot electrons in bulk maiepialst^Tca and image-potential 
states at metal surfacescjEj have been probed. 

These new and powerful experimental techniques, 
based on high resolution direct and inverse photoemis- 
sion as well as time-resolved measurements, have ad- 
dressed aspects related to the lifetime of excited electrons 
and have raised many fundamental questions. The ultra- 
fast laser technology has allowed to probe fast events at 
surfaces in real time and, therefore, extract information 
about elementary electronic processes (with time scales 
from pico to femtoseconds) that are relevant for potential 
technological applications. In general, the two-photon 
photoemission spectroscopy is sensitive to changes of ge- 
ometries, local work functions, and surface potentials 
during layer formation. The interaction of excited elec- 
trons and the underlying substrate governs the cross- 
section and branching ratios of all electronically induced 
adsorbate reactions at surfaces, such as dissociation or 
desorption, and influences the reactivity of the surfaces 
as well as the kinetics of growth. tj Hot-electron life- 
times have long been invoked to give valuable information 
about these processes. 

Inelastic lifetimes of excited electrons with energies 
larger than ~ 1 eV above the Fermi level can be at- 
tributed to electron-electron (e-e) inelastic scattering, 
other processes such as electron-phonon and electron- 
imperfectioii—interactions being, in general, of minor 
importance.Ell A self-consistent calculation of the inter- 



1 



action of low-energy electrons with axuelectron gas was 
first carried out by Quinn and Ferrell.H They performed 
a self-energy calculation of e-e scattering rates near the 
Fermi surface, and derived a formula for the inelastic 
lifetime of hot electrons that is exact in the high-density 
limit. These free-eiectron-gas (FEG) calculations were 
extended by Ritchia2j and QuinrcJ to include, within the 
first-Born and random-phase approximatipns, energijes. 
away from the Fermi surface, and by AdleJlil and Quinnc3 
to take account of the effects of the presence of a periodic 
lattice andj-4p particular, the effect of virtual interband 
transitions.eS Since then, several FEG calculations of e- 
e scattering rates have been performed, rW|itii| inclusion 
of exchange and correlation j-DCC) effectspJEZl chemical 
potential renormaUg|ation,c3EJ plasmon damping,E2l and 
core polarizabilitytil In the case of free-electron materi- 
als, such as aluminum, valence electrons were described 
within the FEG model and atomic generalized oscilku. 
tor strengths were used for inner-shell ionizationOS 
For the description of the IMEP in non-free-electron 
metals, Krolikowski and Spicei£3 employed a semiem- 
pirical approach to calculate the energy dependence of 
the IMFP from the knowledge of density-of-state dis- 
tributions, which had been deduced from photoeiectron 
energy-distribution measurements. Tung et alc^ used 
a statistical approximation, assuming that the inelastic 
scattering of an electron in a given volume clement of the 
solid can be represented by the scattering appropriate to 
a FEG with the electron density in that volume element. 
This approximation was found to predict IMFPs for elec- 
trons in Al that are in good agreement with predictions 
from an electron gas model Ipivs atomic inner-shell con- 
tributions, and these authorf£3 went further to evaluate 
IMFPs and energy losses in various noble and tpLn^ion 
metals. Later on, new methods were proposedE^L^ for 
calculating the IMFP, which were based on a model di- 
electric function whose form was motivated by the use 
of optical data. Though high-energy electpin mean free 
paths now seem to be well understoodj£j~L3 in the low- 
energy domain electrons are more sensitive to the details 
of the band structure of the solid, and a treatment of 
the electron dynamics that fully includes band structure 
effects is necessary for quantitative comparisons with ex- 
perimentally determined attenuation lengths and relax- 
ation times. Ab initio calculations of these quantities in 
which both the electronic Bloch states of the probe elec- 
tron and the dielectric response function of the medium 
are described fror n firs t principles have been performed 
only very recentlyE^CJ 

The self-energy formalism first introduced by Quinn 
and Ferrell for the description of the lifetime of hot 
electrons in a homagatieous electron gas was extended 
by Echenique et aZE3~EZl to quaplitatively evaluate the 
lifetime of image^pptjential statesE3L3 at metal surfaces. 
Echenique et a/E3~Ell used hydrogenic-like states to de- 
scribe the image-state wave functions, they introduced a 
step model potential to calculate the bulk final-state wave 
functions, and used simplified free-electron-gas models to 



approximate the screened Coulomb interactiipn. A three 
band model was used by Gao and LundqvistLJ to describe 
the band structure of the (111) surfaces of copper and 
nickel. They calculated, in terms of Auger transitions, 
the decay of the first image state on these surfaces to the 
n = crystal-induced surface state, neglecting screen- 
ing effects. Self-consistent calculations of the linewidths 
of image-states on copper surfaces have been reported 
recentlyJIO e3 and good agreement with experimentally 
determined decay times has been found. These calcula- 
tions were performed by going beyond a free-electron de- 
scription of the metal surface. Single-particle wave func- 
tions were obtained by solving the Schrodinger ec»ation 
with a realistic one-dimensional model potential,EJ and 
the screened interaction was evaluated in the random- 
phase approximation (RPA). 

This review includes an overview of inelastic lifetimes 
of low-energy electrons in the bulk and at the surface 
of solid materials, as derived within the first-Born ap- 
proximation or, equivalently, linear response theory. In 
the framework of linear response theory, the inelastic en- 
ergy broadening or lifetime-width of probe particles in- 
teracting with matter is found to be proportional to the 
square of the probe charge. Extensions that include the 
quadratic response to externaLppjturbations have been 
discussed by various authors,E3E3 in order to give ac- 
count of the existing dependence of the energii-lass and 
the IMFP on the sign of the projectile charge.EJ'Ej 

Section II is devoted to the study of electron scatter- 
ing processes in a homogeneous electron gas, employing 
various representations of the electronic response of the 
medium. In section III, a general self-energy formulation 
appropriate for the description of inhomogeneous many- 
body systems is introduced. This formulation is applied 
in sections IV and V to review theoretical investigations 
of lifetimes of both hot electrons in bulk materials and 
image-potential states at metal surfaces. Future work is 
addressed in Section VI. 

Unless otherwise is stated, atomic units are used 
throughout, i.e., = h = rUf. = 1. The atomic unit 
of length is the Bohr radius, qq — /m\ = 0.529A, the 
atomic unit of energy is the Hartree, 1 Hartree = e^/ag = 
27.2 eV, and the atomic unit of velocity is the Bohr ve- 
locity, Wo = a c = 2.19 X lO^cm s^^, a and c being the fine 
structure constant and the velocity of light, respectively. 

II. SCATTERING THEORY APPROACH 

We take a homogeneous system of interacting elec- 
trons, and consider an excited electron interacting 
through individual collisions with electrons in the Fermi 
sea. Hence, we calculate the probability P/'J per unit 
time corresponding to the process by which the probe 
particle is scattered from a state 4>i{v) of energy Ei to 
some other state (j^fij) of energy Ef, by carrying one 
electron of the Fermi sea from an initial state 4>i'{r) of 
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energy Ei^ to a final state (f>f'{r) of energy Efi, accord- 
ing to a dynamic screened interaction W {r — r' ; Ei — E f) 
(see Fig. 1). By using the 'golden rule' of time-dependent 
perturbation theory and keeping only tpms of first order 
in the screened interaction, one writesM 



27r 



/,/' 



[W{r-r';E,-Ef)YJ, 



xSiE,-Ef+E'^-E'f), 



(1) 



where 

[W{r- 



dr / dr'0*(r)0*(r')l^(r-r';c.)0/(r)0/,(r')- (2) 



(3) 



Using plane waves for all initial and final states, 
. , > 1 



n 



with energy Wk — and fl being the normalization 

volume, one finds 



(4) 



where Wq^i^ represents the Fourier transform of the 
screened interaction W{r — r';u!). The Kroenecker delta 
and the Dirac delta function on the right-hand side of 
Eq. (Q) allow for wave-vector and energy conservation, 
respectively. 

By summing the probabilities P/ of Eq. (^) over all 
possible states {k\ < qp, qp being the Fermi momen- 
tum), {k'j: > qp) and k/, and noting that each allowed 
k,^ leads to two one-electron states (one for each spin), the 
total scattering rate of the probe electron in the state k^ 
is found to be given by the following expression: 



5]'^|W^q,.Pnk;(l-nk^+q) 



where 



An 

q k; 

XS{UJ ~ CJkJ+q + Ci^kj), 

nk = 0{qF - |k|) 



(5) 



(6) 



represents the occupation number. We have set the en- 
ergy transfer to — LUk_. —uj^.^q, and the prime in the sum- 
mation indicates that the momentum transfer is subject 
to the condition < oj < Wki ^ Ep {Ep is the Fermi 
energy), accounting for the fact that the probe electron 
cannot make transitions to occupied states in the Fermi 
sea. 

With the interaction Wq^ui described by the bare 
Coulomb interaction, that is, Wq^ui — Vq, the summation 



over q in Eq. (g) would be severely divergent, thereby 
resulting in an infinite damping rate. Instead, we assume 
that the Coulomb interaction is dynamically screened, 



(7) 



where ja taken to be the dielectric function of the 
medium .Mej 

For [x)_^_Q, the imaginary part of the RPA dielectric 
functionEIrL^ is simply a measure of the number of states 
available for real transitions involving a given momentum 
transfer q and energy transfer uj: 

Ime^^^ = 27rr!-l Vq ^ 77k(l - 71k+q) 

k 

x5(LJ-CJk+q + t^k)- (8) 

In the limit that the volume of the system Q becomes 
infinite, one can replace sums over states by integrals 
with the following relation 



k 



/(k) 



n 



{2n) 



dk/(k). 



(9) 



and after introduction of Eq. (||) into Eq. (^) , one finds 

(10) 



, f dq Ime^,^^ 



where the prime in the integration indicates that the mo- 
mentum transfer q is subject to the same condition as in 
Eq. (||). With the screened interaction Wq,^; of Eq. (|^) 
described within RPA, one writes 



= 2 



dq 



1 Im [- 



with 

^q,u; - 



. RPA 



(2^)3 'I 

being the RPA dielectric function,0'0 i 



(11) 



In the more general scenario of manx-body theory 
and within the first Born approximation,0 one finds the 
damping rate of an excited electron in the state k^ to also 
be given by Eq. (11), but with the exact inverse dielectric 
function eqjj, as defined in Appendix A. This is the re- 
sult obtaiufid independently by Ouinn and FerrelEa and 
by RitchieLJ. Quinn and FerrelH demonstrated, within 
a self-energy formalism, that the damping rate of holes 
below the Fermi level is also given by Eq. ( pT| ) , with the 
energy transfer uj = Wk^-q — cok- and with the prime in 
the integration indicating that the momentum transfer 
q is subject to the condition < a; < Ep — cjki- For 
small values of jo^k^ — Ep\, holes inside the Fermi sea 
(wki < Ep) are found to damp out in the same way as 
electrons outside (wk; > Ep), as shown in Fig. 2. 

If one is to go beyond RPA and introduce, through the 
factor (1 — Gq.i^) (see Appendix A), the reduction in the 
e-e interaction due to the existence of a local XC hole 
around electrons in the Fermi sea, the dielectric function 
entering Eq. (fll) is 
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1 + 



RPA 



1 



1 



G (f^ 



RPA 



1 



(12) =-2E'/ / dr'^*(r)0}(r') 



where Gq.i^ is the 



lUj-caUed local-field factor, first intro- 
duced by Hubbard .EH 

If one accounts, through the factor (1 — Gq.cj), for the 
existence of a local XC hole around electrons in the Fermi 
sea and also around the probe electron, the dielectric 
function entering Eq. j|S |the so-called test_charge- 

electron dielectric functionESo (see Appendix A): 



RPA 



RPA 



(13) 



Finally, we note that the inelastic mean free path 
(IMFP) is directly connected to the lifetime r through 
the relation 

\ = VT. (14) 



III. SELF-ENERGY FORMALISM 

In the framework of many-body theory,@ the damping 
rate of an electron with energy Si > Ep is obtained from 
the imaginary part of the electron self-energy: 



-2 Jdr Jdr' 0*(r)ImE(r,r';eO'/'.(r'), 



(15) 



where (pi (r) represents a suitably chosen one-electron or- 
bital of energy Si (see Appen dix . B) ,. 

In the GW approximationJ12Jll22l one considers only 
the first-order term in a series expansion of the self- 
energy in terms of the screened interaction W{r,r',uj). 
This is related to the density-response function x(r, r', tf ) 
of Eq. (A2), as follows 



W{r,r';Lu) = v{r — r') + J dri J dr2u(r — ri) 

xx(ri,r2,w)w(r2 - r'), (16) 

where v{r — r') represents the bare Coulomb potential. 

Within RPA, the density-response function satisfies 
and integral equation (see Eq. (A6)), and is obtained 
from the knowledge of the density-response function of 
noninteracting electrons. If, to the same order of ap- 
proximation, one replaces the exact one-particle Green 
function by its noninteracting counterpart, the imaginary 
part of the self-energy can be evaluated explicitly: 

ImS(r,r';£,; > ^;^) = ^ V}(r')ImW^(r,r';c^),/)/(r), 
/ 

(17) 

where w = — £/, and the prime in the summation in- 
dicates that states / (r) available for real transitions are 
subject to the condition that < w < — Ep- Introduc- 
tion of Eq. (|l^ into Eq. (|l|) yields 



xlmW{r,r';uj)Mr'Wf(r). (18) 



In the so-called GWF approximation,E23T£j which in- 
cludes XC effects not present in the GW-RPA, the self- 
energy and damping rate of the excited electron are of 
the GW form, i.e., they are given by Eqs. ( p^ ) and (|l^), 
respectively, but with an effective screened interaction 

W{r,r';u)) — v{r — r') + J dri J dr2 [t;(r — ri) 

+ X^^(r,ri)]x(ri,r2,c.)z;(r2-r'), (19) 

the density-response function now being given by Eq. 
(^). The kernel K'"'{r,r') entering Eqs. (||) and (^) 
accounts for the reduction in the e-e interaction due to 
the existence of short-range XC effects associated to the 
probe electron and to screening electrons, respectively. 

A. Homogeneous electron gas 

In the case of a homogeneous electron gas, single- 
particle wave functions are simply plane waves, as defined 
in Eq. (||). By introducing these orbitals into Eq. ([T^), 
the damping rate of an electron in the state is found to 
be given by Eq. ( 11 ) with the dielectric function of either 
Eq. ( p^ ) or Eq. ( 13 ) , depending on weather the screened 
interaction of Eq. (16) or Eq. (|l9|) is taken in combinaj 
tion with the density-response function of Eq. ( |A8|) .E£! 
This is an expected result, since these calculations have 
all been performed to lowest order in the screened inter- 
action. 



B. Bounded electron gas 

In the case of a bounded electron gas that is trans- 
lationally invariant in the plane of the surface, single- 
particle wave functions are of the form 



'Ak||,z(r) 



VA 



with energies 



^kii ,i 



lf2 

2 



(20) 



(21) 



where the z-axis has been taken to be perpendicular to 
the surface. Hence, the wave functions (f>i{z) and energies 
Si describe motion normal to the surface, k|| is a wave 
vector parallel to the surface, and A is the normalization 
area. 

Introduction of Eq. (|^) into Eqs. (|l|) and (^ yields 
the following expressions for the damping rate of an elec- 
tron in the state </'k||,j(r) with energy ekii.j: 
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-1 - -2 / dz / dz' f 4^0* (z) 



(2^)2 

xI^lE(z,z';q||,ek||,^)0^(2;') (22) 



and 



xIniVK(z, z'; q,,, c^)0/(z)0,(z'), (23) 

respectively, where uj = Eky,! — eky-qn,/- Here, 
i;(2:,z';q||,w) and W^(z, z'; q|| , w) represent the two- 
dimensional Fourier transforms of the electron self-energy 
S(r,r';ti;) and the screened interaction W^(r, r'; w). 

C. Periodic crystals 

For periodic crystals, single-particle wave functions are 
Bloch states 
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4k,*(r) = --=e"'''uu.»(r). 



(24) 



and one may introduce the following Fourier expansion 
of the screened interaction: 

Wir, r'; to) ^ Q-^ f dq ^ ^ e'(i+^) '- 

xe-'('i+°')""VG,G'(q,t^), (25) 

where the integration over q is extended over the first 
Brillouin zone (BZ), and the vectors G and G' are recip- 
rocal lattice vectors. Introducing this Fourier represen- 
tation into Eq. (^8|), one finds the following expression 
for the damping rate of an electron in the state (/)k,i(r) 
with energy Sk^f. 

xB,/(q+G')ImM^G,G'(q,^), (26) 



or, equivalently, 



G G' 



B*.(q+G)i?./(q + G') 



BZ 



|q+G|- 



Within RPA, one writes 

eG,G'(q,w) = ^G,G' - WG(q)XG,G'(q'^)' (29) 

where vc (q) represent the Fourier coefficients of the bare 
Coulomb potential, 



WG(q) 



47r 



|q+G|2' 



(30) 



and Xg g' (q^ ^) Fourier coefficients of the 

density-response function of noninteracting electrons. 



x°G.G'(q,^) = f^-7 ^'^EE 



/k,n /k-t-q,i 



ek,n - £k+q,ri' + (w + i??) 

x(0k,„|e-'(''+°)--|0k+q,„') 

x(0k+q,„'|e'''»+°'^-"|0k,„), (31) 

where 77 is a positive infinitesimal. The sums run over 
the band structure for each wave wave vector k in the 
first BZ, and /k,„ are Fermi factors 



fk.n — 0{Ef — £Vi,n)- 



(32) 



Couplings of the wave vector q -f G to wave vectors 
q -I- G' with G 7^ G' appear as a consequence of the 
existence of electron-density variations in real solids. If 
these terms, representing the so-called crystalline local- 
field effects, are neglected, one can write 



|i?./(q + G)r 



|q+G| 



Ini[eG,G(q,^)" 
|eG,G(q,w)P 



(33) 



xlm 



-eG^G'(q.^' 



(27) 



The imaginary part of eG,G(q, ^) represents a measure of 
the number of states available for real transitions involv- 
ing a given momentum and energy transfer q-l- G and w, 
respectively, and the factor |eG,G(q, w)|~^ accounts for 
the screening in the interaction with the probe electron. 
Initial and final states of the probe electron enter through 
the coefficients i?j/(q + G). 

If one further replaces in Eq. ( |3^ ) the probe electron 
initial and final states by plane waves, and the matrix 
coefficients eG,G(q, w) by the dielectric function of a ho- 
mogeneous electron gas. 



where uj = ek.i — Ek-q./, and 

B^M)^ I dr0*,(r)e'i-0k-q./(r). 



eG,G(q,t^) e(|q-f G|,w), 



(34) 



(28) 



W^G,G'(q, w) are the Fourier coefficients of the screened 
interaction, and eQ^Q,(q, w) are the Fourier coefficients 
of the inverse dielectric function. 



then Eq. (33) yields the damping rate of excited electrons 
in a FEG, as given by Eq. (|ll|). 

We note that the hot-electron decay in real solids de- 
pends on both the wave vector k and the band index 
i of the initial Bloch state. As a result of the symme- 
try of these states, one finds that T~^(S'k, i) = T^^(k, z). 
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with S representing a point group symmetry operation 
in the periodic crystal Hence, for each value of the hot- 
electron energy the scattering rate t~^{E) is defined by 
averaging T~^(k, z) over all wave vectors k lying in the 
irreducible element of the Brillouin zone (IBZ), with the 
same energy, and also over the band structure for each 
wave vector. 



IV. LIFETIMES OF HOT ELECTRONS IN 
METALS 

A. Jellium model 

Early calculations of inelastic lifetimes and mean free 
paths of excited electrons in metals were based on the 
'jellium' model of the solid. Within this model, valence 
electrons are described by a homogeneous assembly of 
electrons immersed in a uniform background of positive 
charge and volume 17. The only parameter in this model 
is the valence-electron density uq, which we represent in 
terms of the so-called electron-density parameter de- 
fined by the relation 1/no = (4/3) vr {rg ao)^, uq being the 
Bohr radius. Hence, the damping rate of a hot electron 
of energy E — uj^. is obtained, within this model, from 
Eq. (O) with lu — ujk.. — uj\^.^q. 

In the high-density limit {vg 0), XC effects as well as 
high-order terms in the expansion of the scattering prob- 
ability in terms of the screened interaction are negligible. 
Thus, in this limit the damping rate of hot electrons is 
obtained from Eq. (11) with use of the RPA dielectric 
function. 

Now we focus on the scattering of hot electrons just 
above the Fermi level, i. e., E — Ep << Ep. As the 
energy transfer lu cannot exceed the value E — Ep , the 
frequency entering Im [cq J^] is always small, one can take 



Im [- 



Iml 



3 ^q.o'^' 



(35) 



and the {E — EpY quadratic scaling of the hot-electron 
damping rate is predicted. If one further replaces, within 
the high-density limit {qp — > oo), the static dielectric 
function e^ Q by the Thomas-Fermi approximation, and 
extends, at the same time, the maximum momentum 
transfer {q ^ 2 qp) to infinity, then one finds 



(TT/gf )3/2 {E - Epf 
16 ki 



(36) 



If we replace ki — > qw-V^ Eq. (p6|), then the damping rate 
of Quinn and Ferrel£3 is obtaine d. T q^- as given by Eq. 
(C8). For the lifetime, one writesll2Zl 



TQF = 263 r^T^/^ i^E ~ Ep)-'^iseY^ 



(37) 



denominator |eq.oP accounts for the screening in the in- 
teraction between the hot electron and the Fermi sea. 
Hence, the hot-electron lifetime is determined by the 
competition between transitions and screening. Though 
increasing the electron density makes the density of states 
(DOS) larger, momentum and energy conservation pre- 
vents, in the case of a FEG, the sum over available states 
from any dependence on Ts, as shown by Eq. (|35|). As a 
result, the scattering rate of hot electrons in a FEG only 
depends on the electron-density parameter through the 
screening and the initial momentum k^. High densities 
make the interaction weaker [the integration of |eq,(j|~^ 
scales, in the high-density limit, as qp'^^'^] and momenta 
of excited electrons larger [1/ki qF^]^ which results in 



the 



-5/2 



scaling described by Eq. (^ 



In Eq. (^) , Im eq represents a measure of the num- 
ber of states available for real transitions, whereas the 



In Fig. 3 we represent the ratio t/tqp, versus E — Ep, 
for an electron density equal to that of valence electrons 
in copper (ts = 2.67), as obtained from Eq. ( pT| ) with 
the full RPA dielectric function (solid line) and from Eq. 
( ^ ) (dashed hue). Though in the hmit E Ep the 
available phase space for real transitions is simply E — 
Ep, which yields the {E — EpY quadratic scaling of Eqs. 
( ^ ) and (^^, as the energy increases momentum and 
energy conservation prevents the available phase space 
from being as large as E — Ep. Hence, the actual lifetime 
departures from the ki/{E — Ep)^ scaling predicted for 
electrons in the vicinity of the Fermi surface, differences 
between full RPA calculations (solid line) and the results 
predicted by Eq. (pfl) (dashed line) ranging from ^ 2% 
at i5 « Ep to ^ 35% at E — Ep = 5 eV. For comparison, 
also represented in this figure is the ratio t/tqp obtained 
from the approximations of Eqs. (C4) (dotted line) and 
(dashed-dotted line). 
The result of goinjg-bci/ond the RPA has been discussed 
by variojis authors .c3l2I In an early paper, Ritchie and 
Ashleyc3 investigated the simplest exchange process in 
the scattering between the probe electron and the elec- 
tron gas. Though this exchange contribution to the e-e 
scattering rate is of a higher order in the electron-density 
parameter than the direct term, it was found to yield, 
for Ts = 2.07 and E ~ Ep, a ~ 70% increase with re- 
spect to the RPA lifetime, and an even larger increase in 
the case of metals with Vg > 2. This reduction of the e-e 
scattering rate appears as a consequence of the exclusion 
principle keeping two electrons of parallel spin away from 
the same point, thereby reducing their effective interac- 
tion. 

Neither the effect of Coulomb correlations between the 
probe electron and the electron gas, which also influence 
the e-e mutual interaction, nor XC effects between pairs 
of electrons jidthin the Femii sea were included by Ritchie 
and Ashley.c3 Kleinmano included not only XC between 
the incoming electron and an electron from the Fermi sea 
but also XC between pairs of electrons within the Fermi 
sea, and found a result which reduced the ~ 70% increase 
obtained by Ritchie and Ashley for Al to a ~ 1% in- 
crease. Alternative approximations for the XC corrected 
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e-e interaction Wjeije derived by Penned and by Kukkonen 
and Overhauser.Cj From an evaluation of the test_charge- 
electron dielectric function of Eq. ( p^ ) and with use of 
a static local- field factor, Penncll concluded that the in- 
troduction of exchange and correlation has little effect 
on the lifetime of hot elecjttpns, in agreement with early 
calculations by Kleinman.E3 

As we are interested in the low-frequency (w — > 0) be- 
haviour of the electron gas, we can safely approximate 
the local- field factor by t he sta tic limit, Gq^o, which we 
choose to be given by Eg. ( Al5| ). Our results, as obtained 
from Eq. (ph with the dielectric function of either Eq. 
( p^ or Eq. (|13D are presented in Figs. 4 and 5 by dashed 
and dotted lines, respectively, as a function of rg for hot 
electrons with E — Ep ~ 1 eV (Fig. 4), and as a function 
of E — Ep with 7's — 2.67 (Fig. 5). Solid lines represent 
RPA calculations, as obtained with the local-field factor 
Gq,ui set equal to zero. We note from these figures that 
local-field corrections in the screening reduce the lifetime 
of hot electrons in a FEG with an electron density equal 
to that of valence (4s^) electrons in Cu (rg = 2.67) by 
~ 20%. However, this reduction is slightly more than 
compensated by the large enhancement of the lifetime 
produced by the existence of local-field corrections in the 
interaction between the probe electron and the electron 
gas. As a consequence, RPA calculations (solid line) pro- 
duce lifetimes that are shorter than more realistic results 
obtained with full inclusion of XC effects (dotted line) by 
~ 5%. 

Instead of calculating the damping rate r^^ on-the- 
energy-shell {E = uj\^. ) , Lundqvistca expanded the elec- 
tron self-energy in the deviation of the actual excitation 
energy E from the independent-particle result, showing 
that near the energy-shell {E ~ WkJ interactions renor- 
malize the damping rate by the so-called renormaliza- 
tion constant Z\^.. Based on Lundqvist's calculations, 
SheltonO derived IMFPs for various values of and 
for electrons with energies between Ep and ^ 2b Ep. 
The resulting IMFPs were larger than those obtained by 
QuinnO by roughly 5 — 20%, depending on and the 
electron energy. 

In the case of excited electrons near the Fermi level the 
renormalization constant, as obtained within the GW- 
RPA, is nearly real and fc-independent. In the metallic 
density range (r^ ~ 2 — 6) one finds Z ~ 0.8 — 0.7, and 
the resulting lifetimes are, therefore, larger than those 
obtained from Eq. (|l]) by - 20%. 



B. Statistical approximations 

In order to account for the inelastiCpScattering rates of 
non-free electron materials, Tung et a/Eil applied a statis^ 
tical approximation first developed by Lindhard et al .E£j 
This approximation is based on the assumption that the 
inelastic electron scattering of electrons in a small vol- 
ume element dr at r is the same as that of electrons in a 



FEG with density equal to the local density. 

Within the statistical approximation of Ref. for a 
given density distribution n(r) one finds the total scatter- 
ing rate by averaging the corresponding local quan- 
tity [^(r)] over the volume of the solid: 



dr T [n(r)] , 



(38) 



By calculating spherically symmetric -electron density 
distributions n(r) in a Wigner-Seitz celljl£j the total scat- 
tering rate is obtained from 



1 

ws 



drr^T ^ 



(39) 



where flws and ryys represent the volume and the radius 
of the Wigner-Seitz sphere of the solid. 

Alternatively, foUowijig the idea of using optical data 
in IMFP |-GalGulationsJj a number of approaches were 
developedE3~E3 to compute a model energy-loss function 
Im [— EqL] foi' r6al solids and then obtain inelastic scat- 
tering rates from Eq. (11). In these approaches the 
model energy-loss function is set in the limit of zero wave 
vector equal to the imaginary part of the measured op- 
tical inverse dielectric function,Eij Im and it is 
then extended into the non-zero wave vector region by a 
physically motivated recipe. 



Combining the statistical method of Ref. 54 with the 
use of optical data, PenrO developed an improved al- 
gorithm to evaluate the dielectric function of the mate- 
rial. The Penn algorithm is based on a model dielectric 
function in which the momentum dependence is deter- 
mined by averaging the energy-loss function of a FEG, 
Im [-l/eq^f^] , as follows 

nOO 

Im [-eq = / diop G(c.p)Im [l / e^^f (cOp)] , (40) 
Jo 



where 



G(c.) = — Im[-l/eJ] 

TTLU 



(41) 



[Xhe Penn algorithm has been employed by Tanuma et 
aM to calculate IMFPs for 50 to 2000 eV in a variety 
of materials comprising elements, inorganic compounds, 
and organic compounds. Recently, several other groupsEll 
have calculated IMFPs from optical data in a manner 
similar to that proposed by PennJ£3 with some differences 
in approach^nd high-energy IMFPs now seem to be well 
understood^. 

The effect of d-electrons in noble anetals has been re- 
cently investigated by Zarate et al,cl}i by including the 
d-band contribution to the measured optical inverse di- 
electric function into a FEG description of the s — p part 
of the response. 

In order to account for the actual DOS in real mate- 
rials, early IMFP calculations were carried out by Kro- 
likowski and Spicer,E2l with the explicit assumption that 



7 



the matrix elements of the screened e-e interaction enter- 
ing Eq. d^) are momentuiB-4«jdependent. This so-called 
"random" k approximatioiilli3-ti3 has proved to be useful 



in cases where the DOS plays a key role in the determi- 
nation of scattering rates, as in the case of ferromagnetic 
materials, Eialli^ thereby allowing to explaiati|e existence 
of spin-dependent hot-electron lifetimesJHa'EZi Although 
this method, due to its simplicity, cannot provide full 
quantitative agreement with the experiment, it provides 
a useful tool for the analysis of experimental data, thus 
allowing to isolate the effects that are directly related to 
the DOS. 

Figs. 6 and 7 show the lifetime versus energy, for repre- 
sentative free-electron-like and non-free-electron-like ma- 
terials, Al (Fig. 6) and Cu (Fig. 7), respectively. First of 
all, we consider a relatively free-electron-like solid such 
as Al (see Fig. 6). The contribution to the inelastic scat- 
tering of low-energy electrons in Al coming from the ex- 
citation of core electrons is negligible. Hence, statistical 
approximations yield results that nearly coincide with the 
FEG calculation with rs = 2.07. However, the effective 
number of valence electrons in Al is 3.1 rather than 3 (the 
actual number of valence electrons), and lifetimes calcu- 
lated from the statistical model of Ref. 54 are, therefore. 



the model energy-loss function of Eq. ( |40| ) with either 
the recipe described by Salvat et al (dashed line) or with 
the measured optical response function taken from Ref. 



lie (dashed-dotted-dotted-dotted line), and the empiri- 



slightly larger than those obtained within a FEG descrip- 
tion. At higher energies, new contributions to the inelas- 
tic scattering come from the excitation of core electrons, 
and FEG lifetimes would, therefore, be much longer than 
those obtained from the more realistic statistical approx- 
imations. 

For non-free-electron-like materials such as Cu, the 
role of d states in the electron relaxation process is of 
crucial importance, even in the case of very-low-energy 
electrons. The effective number of valence electrons in 
Cu that contribute through the average of Eq. (^, at 
low electron energies, to the inelastic scattering ranges 
from 2.5 far from atomic positions to ~ 7.5 in a region 
where the binding energy is already too large. Since an 
enhanced electron density results in a stronger screening 
and, therefore, a longer lifetime (see, e.g., Eq. (37)), the 
statistical approximation yields lifetimes that are longer 
than those obtained within a FEG model with the elec- 
tron density equal to that of valence (4s^) electrons in 
Cu {rg = 2.67), but shorter than those obtained within a 
FEG model with the electron density equal to that of all 
4s^ aud 3(i^° electrons in Cu. We note that the theory of 
PcpaEj gives shorter lifetimes than the theory of Tung et 
al j£3 which is the result of spurious contributions to the 
average energy- loss function of Eq. ( ^0| ) from Im [l/ej^*'*] 
at very low-frequencies. 

For comparison with the 'universal' relationship t^^ = 
0.13 (i? — Ep) proposed by Goldmann et al[3 for Cu, 
on the basis of experimental angle-resolved inverse pho- 
toemission spectra, lifetime- widths of high-energy 
electrons in Cu are represented in Fig. 8. Solid and 
dashed-dotted lines represent results obtained from Eqs. 
( pi] ) and (|3^). Dashed and dashed-dotted-dotted-dotted 
lines represent the result of introducing into Eq. (|l] 



cal formula of Goldmann et alt3 is represented by a dot- 
ted line. We note from this figure that while at low elec- 
tron energies increases quadratically with E — Ep, 
a combination of inner-shell and plasmon contributions 
results in lifetime-widths that approximately reproduce, 
for electron energies in the range ~ LO.— 50 eV above the 
Fermi level, the empirical predictiontj that the lifetime- 
width increases linearly with increasing distance from 
Ep. 

High-energy lifetime-widths and IMFPs seem to be 
well described by model dielectric functions, by assuming 
that the probe wave functions are simply plane waves. 
Nevertheless, in the case of low-energy electrons band 
structure effects are found to be important, even in the 
case of free-electron-like metals such as Al, and a a treat- 
ment of the electron dynamics that fully includes band 
structure effects is necessary for quantitative comparisons 
with the experiment. 



C. First- principles calculations 

Ah initio calculations of the inelastic lifetime of 
hot electronS|-4p-| metals have ibeea carried out only 
very recently.L3E3 In this workJ^'Ej Bloch states were 
first expanded in a plane- wave basis, and the i^Jahft 



Sham equation of density-functional theory (DFT)IH3'EiI 
was then solved by invoking the local-density approx- 
imation (LDA). The electron- ion interaction was de- 
scribed by means-of a non-local, norm-conserving ionic 
pseudopotential,tia and the one-electron Bloch states 
were then used to evaluate both the Bif coefhcients and 
the dielectric matrix ecc entering Eq. (p7|). 

First-principles calculations oL the average lifetime 
t{E) of hot electrons in real A1,L^ as obtained from Eq. 
( p7| ) with full inclusion of crystalline local field effects, 
are presented in Fig. 9 by solid circles. As an Al crys- 
tal does not present strong electron-density gradients 
nor special electron-density directions (bondings), con- 
tributions from the so-called crystalline local-field effects 
are found to be negligible. On the other hand, band- 
structure effects on the imaginary part of the inverse di- 
electric matrix are approximately well described with the 
use of a statistical approximation, as obtained from Eq. 
( ^ ) (dotted line), thereby resulting in hfetimes that are 
just slightly larger than those of hot electrons in a FEG 
with Ts = 2.07 (solid line). Therefore, differences be- 
tween full ah initio calculations (solid circles) and FEG 
calculations (solid line) are mainly due to the sensitiv- 
ity of hot-electron initial and final wave functions on the 
band structure of the crystal. When the hot-electron 
energy is well above the Fermi level, these orbitals are 
very nearly plane-wave states and the lifetime is well de- 
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scribed by FEG calculations. However, in the case of 
hot-electron energies near the Fermi level, initial and fi- 
nal states strongly depend on the actual band structure 
of the crystal. Due to the opening, at these energies, of 
interband transitions, band structure effects tend to de- 
crease the inelastic lifetime by a factor that varies from 
~ 0.65 near the Fermi level {E — Ep = 1 eV) to a factor 
of ~ 0.75 tor E-Ef = 3eV. 

Ab initio calculations-of the average lifetime t{E) of 
hot electrons in real Cu,Eil the most widely studied metal 
by TR-2PPE, are exhibited in Fig. 10 by sohd circles, 
as obtained from Eq. (^7|) with full inclusion of crys- 
talline local field effects and by keeping all 4s^ and Sd^" 
Bloch states as valence electrons in the pseudopotential 
generation. The lifetime of hot electrons in a FEG with 
the electron density equal to that of valence (4s^) elec- 
trons in Cu {vs — 2.67) is represented by a solid line, 
and the statistically averaged lifetime, as obtained from 
Eq. (|39|), is represented by a dotted line. These calcula- 
tions indicate that the lifetime of hot electrons in real Cu 
is, within RPA, larger than that of electrons in a FEG 
with Ts = 2.67, this enhancement varying from a factor 
of ~ 2.5 near the Fermi level {E — Ep — 1.0 eV) to a 
factor of ~ 1.5 for E — Ep = 3.5 eV. Ab initio calcu- 
lations of the lifetime of hot electrons in Cu, obtained 
by just keeping the 4s^ Bloch states as valence electrons 
in the pseudopotential generation, were also performed, 
and they were found to nearly coincide with the FEG 
calculations. Hence, d-band states play a key role in the 
hot-electron decay mechanism. 

In order to address the various aspects of the role 
that localized d-bands play on the lifetime of hot elec- 
trons in Cu, now we neglect crystalline local-field effects 
and present the result of evaluating hot-electron lifetimes 
from Eq. (p^). First, we replace hot-electron initial and 
final states in |_Bi/(q, G)p by plane waves, and the di- 
electric function in |eG,G(q, by that of a FEG with 
rs = 2.67. If we further replaced Im [eG,G(q, w)] by that 
of a FEG, i.e., Eq. (||), then we would obtain the FEG 
calculation represented by a solid line. Instead, Ogawa et 
alt3 included the effect of d-bands on the lifetime by com- 
puting the actual number of states available for real tran- 
sitions, within Eq. (^), from the band structure of Cu, 
and they obtained a result ti^t is for E ~ Ep > 2 eV well 
below the FEG calculation.ti3 However, if one takes into 
account, within a full description of the band structure of 
the crystal in the evaluation of Im [eG,G(q, i^)] (see Eqs. 
( p9| ) and (|3l|)), couplings between the states participating 
in real transitions, then one obtains the result represented 
in Fig. 10 by open circles. Since the states just below 
the Fermi level, which are available for real transitions, 
are not those of free-electron states, localization results 
in lifetimes of hot electrons with E — Ep < 2 eV (open 
circles) that are slightly larger than predicted within the 
FEG model of the metal. At larger energies this band- 
structure calculation (open circles) predicts a lower life- 
time than within the FEG model, due to opening of the 



d-band scattering channel dominating the DOS with en- 
ergies from ~ 2 eV below the Fermi level. Thus, this cal- 
culation shows at E — Ep ^ 2 eV a slight deviation from 
the quadratic scaling predicted within the FEG model, 
in qualitative agreement with experimentally determined 
decay times in Cu. 

While the excitation of d electrons diminishes the life- 
time of hot electrons with energies E — Ep > 2 eV, d 
electrons also give rise to additional screening, thus in- 
creasing the lifetime of all hot electrons above the Fermi 
level. That this is the case is obvious from the band- 
structure calculation exhibited by full triangles in Fig. 
10. This calculation is the result obtained from Eq. ( |3^ ) 
by still replacing hot-electron initial and final states in 
|i?o/(q + G)|^ by plane waves (plane- wave calculation) 
but including the full band structure of the crystal in the 
evaluation of both Im [eG,G(q, t^)] and |eG,G(q, ^^)|~^- 
The effect of virtual interband transitions giving rise to 
additional screening is to increase, for hot-electron en- 
ergies under study, the lifetime by a factor of « 3, in 
qualitative agreement with the approximate prediction 
of Quinnc3 and with the use of the statistical average of 
Ref. ||. 

Finally, band-structure effects on hot-electron energies 
and wave functions are investigated. Full band-structure 
calculations of Eq. (|27| ) with and without (see also Eq. 
(^3|)) the inclusion of crystalline local field corrections 
were carried out,tJ and these corrections were found to 
be negligible for E — Ep > 1.5 eV, while for energies very 
near the Fermi level neglection of these corrections re- 
sulted in an overestimation of the lifetime of less than 5%. 
Therefore, differences between the full (solid circles) and 
plane-wave (solid triangles) band-structure calculations 
come from the sensitivity of hot-electron initial and final 
states on the band structure of the crystal. When the 
hot-electron energy is well above the Fermi level, these 
states are very nearly plane-wave states for most of the 
orientations of the wave vector, and the lifetime is well 
described by plane-wave calculations (solid circles and 
triangles nearly coincide ior E — Ep > 2.5 eV). However, 
in the case of hot-electron energies near the Fermi level 
initial and final states strongly depend on the orienta- 
tion of the wave vector and on the shape of the Fermi 
surface. For most orientations, flattening of the Fermi 
surface tends to increase the hot-electron decay rate,cll 
while the existence of the so-called necks on the Fermi 
surface of noble metals results in very small scattering 
rates for a few orientations of the wave vector. After 
averaging r~^(k, n) over all orientations, Fermi surface 
shape effects tend to decrease the inelastic lifetime. 

Scaled lifetimes, t x {E — Ep)'^, of hot electrons in 
Cu are represented in Fig. 11, as a function oi E — Ep. 
Results obtained, within RPA, from Eqs. (|ll|) and ( |39| ) 
are represented by solid and dashed-dotted lines, respec- 
tively, the ab initio calculations of Ref. ^ are repre- 
sented by solid circles, and the dashed line represents 
the calculations described in Ref. 111. These model 
calculationgliJ show that above the d-band threshold, at 
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'-^ — 2eV relative to the Fermi level, d-band electrons 
can only participate in the screening, thereby producing 
longer lifetimes, while at larger energies lower lifetimes 
are expected, due to opening of the c?-band scattering 
channel that dominates the DOS with energies ~ 2eV 
below the Fermi level. For comparison, the empirical 
formula proposed by Goldmann et al is represented by a 
dotted line. 

Scaled lifetimes of hot electmns-in Cu, determined 
from a variety of experiments,LZrE2l are represented in 
Fig. 12, as a function oiE — Ep. The energy dependence 
of both ah initio calculations and experimentally deter- 
mined relaxation times shows deviations with respect to 
the {E — Ep)^"^ scaling. On the other hand, though 
there are large discrepancies among results obtained in 
different laboratories, most experiments give lifetimes 
that are considerably longer than predicted within a free- 
electron description of the metal, in agreement with first- 
principles calculations. Measurements of hot-electron 
lifetimes have alstt-Jjeen. performed for othcc, noble and 



netic solids 
ble| 



transition ipetaJsJl^'tHEj Ej simple metals,E£lp£prromag- 



O and high-Tc superconductors.E^ (see Ta- 



V. LIFETIMES OF IMAGE-POTENTIAL STATES 
AT METAL SURFACES 

A. Concept and development of image states 

A metal surface generates electron states that do not 
exist in a bulk metal. These states can be classified 
into two groups, according to their charge density lo- 
calization relative to the surface atomic layer: intrin- 
sic surface states and image-potential states. The so- 
called intjdxi^ic surface states, predicted by TammEH2l and 

at the surface atomic 
^ appear in the vacuum 



Shockley, 
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are localized m, 
layer. Image-potential state; 
region of metal surfaces with a band gap near the vac- 
uum level, as a result of the self-interaction of the electron 
with the polarization charge it induces at the surface. Far 
from the surface, into the vacuum, this potential well ap- 
proaches the long-range classical image potential, — l/4z, 
z being the distance from the surface, and it gives rise 
to a series of image-potential states localized outside the 
metal. 

In a hydrogenic model, with an infinitely high repulsive 
surface barriet. these states form a Rydberg-like series 



with energieiSI 



En — 



-0.85ey 



(42) 



converging towards the vacuum level E^ — Q. The corre- 
sponding eigenfunctions are given by the radial solutions 
of an s-like state of the hydrogen atom 



</)„(z)cxzi?5r°(^/4) 



(43) 



The lifetime of these states scales 
the quantum number n, as foUowa^ 



Tn OC n 



■asymptotically with 



(44) 



For a finite repulsive surface barrier, as is the case for 
real metal surfaces, Eq. (42) may be transformed into 



Er, 



-OMeV 
(n + a)2 ' 



(45) 



where a is a quantum defect depending on both the 
energy-gap position and width and also, on the position 
of the image state relative to the gap.Eao 

After demonstration of the resplubility of the image- 
state series on metal— ,siiffaces,Ej these states were 



Binding energies of these 
states _iiaYe-|been measured by inverse photoemissioa. 
(IPE),E3"Ea two-photon photoemission (2PPE),Ea"E2l 
and ti|mctj^solved two-photon photoemission (TR- 
2PPE).LJ"E3 These measurements have provided highly 
accurate data of image-state binding energies at the sur- 
faces of many noble and transition metals, as shown, e.g., 
in Ref. Along with the measurements of image-state 
energies, the dispersion of these states has also been mea- 
sured, and it has been found that only on a few surfaces 
such as Ag(lOO), Ag(lll), and Ni(lll) the first image 
state is characterized Jay an effective mass that exceeds 
the free-electron mass.LJ At the same time, theoretical ef- 
forts have been directed to create relatively simple mod- 
els that reproduce the experimentally observed binding 
energies and effective masses of ima^ states, _aiid also 

First- 



to evaluate the image-plane position.^ 
principL 
ried out, I 



13C- 



13i 
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ions of image states have also been car- 
with various degrees of sophistication. 
This intensive work on image states has resulted in an 
understanding of some of the key points of the physics of 
these states and of the relatively extensive data-base of 
their energies on noble and transition metals. 



B. Lifetimes of image states 

1. Introduction 

In contrast to the relatively simple spectroscopic 
problem of determining the position of spectral fea- 
tures that directly reflect the density of states and 
which may be, in principle, calculated within a one- 

pr line- 
These 



electron theory, the study of spectral widtl 



14£ 



shapes is essentially a many-body problem 
spectral widths appear as a result of electron-electron, 
electron-defect, j-electron-phonon, and electron-photon 
interactions,Ed£rii£j and theK|_axe also influenced by 
phonon-phonon interactionsO'E^il Accurate and system- 
atic measurements of the linewidth of image states on 
metal surfaces were carried out with the use of 2PPE 
spectroscopy (for a review see, e.g., Ref. These ex- 
periments gave smaller values for the image-state lifetime 
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than the ones obtained ia recent very-high resolution 
TR-2PPE measurements-CaEl The reason for this dis- 
crepancy is that the 2PPE linewidth contains not only 
an energy relaxation contribution (intrinsic lifetime), 
tributions tha arise from phase-relaxation 



but also cm 
processes, 1^^° 



The first estimatiorO of the lifetime of image states 
used simple wave-function arguments, to show that the 
lifetime of image states asymptotically increases with the 
quantum number n, as in Eq. (|4^). Nearly twenty years 
later, this prediction was confirmed experimentally for 
the (100) surface of Cu, for which lifetimes of the first 
six image states wem measured with the use of quantum- 
beat spectroscopy^ The first quantitative evaluation of 
the lifetime of image states, as obtained within the self- 
energy formalism, was reported in Ref. In this calcu- 
lation hydrogenic-like states were used, with no penetra- 
tion into the solid, to describe the image-state wave func- 
tions, a step model potential was introduced to calculate 
the bulk final-state wave functions, and a simplified free- 
electron-gas (FEG) model was utilized to approximate 
the screened Coulomb interaction. More realistic wave 
functions, allowing for penetration of the electron intallie. 
crystal, were introduced in subsequent calculations Ej'EII 
In these evaluations the linewidth of the first image state 
at the r point was shown to be directly proportional to 
the penetration, and the prediction of Eq. (E3) was con- 
firmed. 

The penetration of an image state into the bulk is de- 
fined as 



bulk 



(46) 



thereby giving a measure of the coupling of this state 
to bulk electronic states. This coupling, weighted by 
the screened interaction, is responsible for the decay of 
image states through electron-hole pair creation. Intu- 
itively, it seems clear that the larger the penetration the 
stronger the coupling and, therefore, the smaller the life- 
time. This idea was exploited to qualitatively lExplain the 
linewidth of image states on various surfaces Jlj and also 
the temperature-dependence of the linewidth of the n=l 
state on Cu(lll).Eij In this heuristic approximation, the 
linewidth of an image state is determined by 



(47) 



where Tb{En) is the linewidth of a bulk state correspond- 
ing to the energy En- The rb{En) value can be obtained 
either from first-principles calculations or from the ex- 
periment. In many angle-resolved photocmission experi- 
ments a linear dependence of the linewidth of bulk s — p 
and d states is observed ior energies in the range 5-50 eV 
above the Fermi leveLEfO 



n{En)^b{Er,~EF) 



(48) 



while the linewidth of bulk states in a FEG shows a 
{En — EfY quadratic scaling for energies near the Fermi 



level, as discussed in section IV. A and in Appendix C. 
Image states on noble and transition metal surfaces have 
energies in the r ang e 4-5 eV above the Fermi level, so that 

for an 
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Eqs. (^) and (|4^) have been applied in Ref. 
estimate of the lifetime broadening, with use of the ex- 
perimentally determined coefficient b = 0.13 for Cu and 
Ag£3 For Au one also uses 0.13, and for Ni and Fe b 
is taken to be O.ISEI and 0.6jiJ respectively. 

Recent TR-2PPE measurements have shown that the 
intrinsic linessddths of the n=l image state on Cu(111)Lj 
and Cu(100)E3't3 are 30 meV and 16.5 meV,r*€spectively, 
while accurate model potential calculationsEll yield pen- 
etrations pi = 0.22 and pi — 0.05, respectively. Ac- 
cordingly, image-state linewidths cannot be explained by 
simply applying Eq. (^7|); instead, contributions to the 
image-state decaying mechanism coming from either the 
evanescent tails of bulk states outside the crystal or the 
existence of intrinsic surface states must also be taken 
into account, together with an accurate description of 
surface screening effects. Here we give the results ob- 
tained within a theory that incorporates these effects 
and that has been used recently to evaluate intrinsic 
linewidths or, eauiyalently, lifetimes of image states on 
metal surfaces. lD~e3 



2. Model potential 



It is well knownEBOlHj that image-state wave func- 
tions lie mainly in the vacuum side of the metal surface, 
the electron moving, therefore, in a region with little po- 
tential variation parallel to the surface. Hence, these 
wave functions can be described, with a reasonable ac- 
curacy, by using a one-dimensional potential that repro- 
duces the key properties of image states, namely, the po- 
sition and width of the energy gap and, also, the binding 
energies of both intrinsic and image-potential states at 
the r point. Such a one-dimensional potential has re- 
cently been proposed for a periodic-film niOjdel with large 
vacuum intervals between the solid films:Eil 



' + Ai cos(27rz/as). 



V{z) = 



z <D 



A20 + A2 cos [I3{z ~ D)], D < z < zi 

A3 exp[-a{z ~ zi)] , zi < z < Zim ^^^^ 

exp [-X{z- Zim)] - 1 



^{z Zijji) 



where the 2;-axis is taken to be perpendicular to the sur- 
face. D is the halfwidth of the film, as is the interlayer 
spacing, Zim represents the image-plane position, and the 
origin is chosen in the middle of the film. This one- 
dimensional potential has ten parameters, Aio, Ai, j42o, 
A2, A3, a, /3, zi. A, and Zim, but only four of them are 
independent, ^lo, Ai, A2, and /3 are chosen as adjustable 
parameters, the other six parameters being determined 



11 



from the requirement of continuity of the potential and 
its first derivative everywhere in space. The parameters 
Ai and reproduce the width and position of the en- 
ergy gap, while A2 and /3 reproduce experimental or first- 
principles energies Eq and Ei of the n — s-p like surface 
state at the T point and the n — I image state, respec- 
tively. This potential is shown schematically in Fig. 
To illustrate the good quality of the image-state wave 
functions obtained with this model potential, we compare 
such wave functions with those obtained with the use 
of first-principles calculations. Probability amplitudes of 
the n = 1 image state on Li(llO), as obtained from ei- 
ther a self-consistent pseudopotential calculation or the 
one-dimensional model potential of Eq. (^) are repre- 
sented in Fig. 14, showing that the agreement between 
the two curves is excellent. The probability amplitude of 
the n = 1 image state on Cu(lOO), as obtained from the 
one-dimensional model potential of Eq. (p9|), also shows 
very good agreement with tbe- result obtained with the 



use of a FLAPW calculationE±3 (see Fig. pj 



Assuming that corrugation effects, i.e., effects asso- 
ciated with spatial variations of the potential in the 
plane parallel to the surface, are not important and that 
the three-dimensional potential can be described by the 
(x, ?/)-plane average, one-electron wave functions and en- 
ergies are taken to be given by Eqs. (^) and (pT]), re- 
spectively. Within a many-body self-energy formalism, 
the linewidth of the n image state with energy ek,™ is 
then obtained from Eq. ( p2| ) or, within either the GW or 
the GWr approximation, from Eq. (ESI). 



3. Results and discussion 

First of all, we present results obtained with use of the 
one-dimensional model potential (MP) of Eq. (^9|), and 
compare with experimental and other theoretical results. 
A summary of experimental results for image-state life- 
times in noble and transition metal surfaces, as obtained 
from either 2PPE or TR-2PPE measurements, is pre- 
sented in Table ||, together with the result of theoretical 
calculations at the F point (k|| = 0). We note that there 
are large differences between 2PPE and TR-2PPE exper- 
imental results for copper and silver surfaces, the lifetime 
broadening derived from recent very-high resolution TR- 
2PPE measurements being smaller than that obtained 
from 2PPE experiments by nearly a factor of 2. 

Theoretical calculations presented in Table |^ can be 
classified into two groups. First, there is the heuristic ap- 
proximation of Eqs. ( [47I ) aud (|48|), which was carried out 
by Fauster and SteinmannO for a variety of metal sur- 
faces. This approach results in a semiquantitative agree- 
ment with 2PPE measurements, except for the (111) sur- 
faces of noble metals and the (100) surface of Ni. Similar 
computations were performed in Ref. ^ for the n — I and 
n = 2 image states on Cu and Ag surfaces, with use of the 
penetration of the image-state wave function that results 



from the one-dimensional model potential of Eq. (|4£ 
Though an accurate description of the penetration of the 
n = 1 image-state wave function yields better agreement, 
in the case of Cu(lll), with the experiment, this heuris- 
tic approach is still in semiquantitative agreement with 
TR-2PPE measurements. A quantitative agreement was 
found for the n — 2 image-state linewidth in Cu(lll). 

In the other group of calculations the many-body self- 
energy formalism described in Section III was-usad for 
the evaluation of the lifetime of image statesJIjiEj re- 
sulting in a quantitative agreement with TR-2PPE mea- 
surements of the lifetime of image states on Cu surfaces 
and showing, therefore, that the present state of the the- 
ory enables to accurately predict the broadening of image 
states on metal surfaces. 

To illustrate the importance that an accurate descrip- 
tion of the self-energy might have on the evaluation of the 
linewidth, we show in Fig. |l| Im [-S(z, z'; k|| = 0, Ei)] 
of the n = 1 image-state electron at the F point 
(k|| = 0) on the (111) and (100) surfaces of Cu. 
Im [— I](z, z'; k|| — 0, -Bi)] is represented in this figure as 
a function of z for a fixed value of z' in the vacuum side of 
the surface (upper panel), within the bulk (middle panel), 
and at the surface (lower panel). It is obvious from this 
figure that-tjhe imaginary part of the self-energy is highly 
nonlocaljl^ and strongly depends on the z and z' coordi- 
nates. We note that Im [— E(z, 2:'; fcy = 0,i?i)] presents 
a maximum at z = z' when z is located at the surface, 
and surface states can, therefore, play an important role 
in the decay mechanism of image states. The magnitude 
of this maximum is plotted, as a function of z', in Fig. 
p^ , showmg that it is an oscillating function of z within 
and reaches its highest value at the surface. 



the bulk 



It is interesting to note from Fig. IGk that the magni 



tude of Im [-S(z, z'; ky = 0, Ei)] is larger at the surface 
for Cu(lll) than for Cu(lOO). Though the (100) sur- 
face of Cu only presents an intrinsic surface resonance, 
in the case of the (111) surface of Cu there is an intrinsic 
surface state just below the Fermi level. This intrinsic 
surface state provides a new channel for the decay of im- 
age states, thereby enhancing the imaginary part of the 
self-energy and the linewidth. The role that the intrinsic 
surface state on Cu(lll) plays in the decaying mecha- 
nism of image states is obvious from Fig. pqb, where 
contributions to the maximum self-energy coming from 
transitions to the intrinsic surface state and from tran- 
sitions to bulk states have been plotted separately by 
dashed and dotted lines, respectively. The intrinsic sur- 
face state provides a ~ 75% of the decay mechanism at 
the surface. The intrinsic-surface-state contribution to 
the total linewidth of the n — IpifM-ge state on Cu(lll) 
was found to be of about 40%.l3l3 Similarly, lower ly- 
ing image states can give noticeable contributions to the 
linewidth of excited, i.e., n = 2,3,... image states. For 
example, the decay from the n = 2 to the n = 1 image 
state on Cu(lOO) yields a linewidth of 0.5 meV, i.e. a 10% 
of the total n = 2 image-state linewidth. The decay from 
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the n = 3 to the n = 1 image state on Cu(lOO) yields a 
linewidth of 0.17 meV, and the decay from the n = 3 to 
the n = 2 image state on Cu(lOO) yields a hnewidth of 
0.05 meV, i.e. ~ 10 - 15% of the total hnewidth. 

Coupling of image states with the crystal occurs 
through the penetration of the image state wave func- 
tion and, also, through the evanescent tails of bulk and 
surface states outside the crystal. To illustrate this point, 
the linewidth F = r^^ can be split as follows 



r = r 



bulk 



(50) 



where Tbuik, ^vac, and Tinter represent bulk, vacuum and 
interface contributions, respectively. They are obtained 
by confining the integrals in Eq. ( ^3| ) to either bulk 
{z < 0, z' < 0), vacuum (z > 0, z' > 0), or vacuum-bulk 
(z < 0, z' > 0; z > 0, z' < 0) coordinates. These separate 
contributions to the line widt h of image states on Cu sur- 
faces are shown in Table III . We note that the contribu- 
tion to r coming from the interference term Tinter is com- 
parable in magnitude and opposite in sign to both bulk 
and vacuum contributions. This is a consequence of the 
behaviour of the imaginary part of the two-dimensional 
Fourier transform of the self-energy, as discussed in Ref. 
|78| . The contributions T^ac and Tinter almost completely 
compensate each other, and, in a first approximation, the 
total linewidth T can be represented by the bulk contri- 
bution Tbuik within an accuracy of ~ 30%. 

The linewidth of image states can vary as a function of 
the two-dimensional momentum k|| . In Fig. ^ we show 
schematically the projection of the bulk band structure 
onto the (111) and (100) surfaces of Cu. In the case of the 
(111) surface of Cu, the n = 1 image state becomes a res- 
onance at |k|| I ^ 0.114 ap^, thereby the image-state wave 
function presenting a larger penetration into the bulk and 
an enhanced linewidth. Table IV shows the linewidth of 
the n = 1 image state on Cu(lll), as obtained foij-tiiree 



different values of ky in the range - 0.114aoMi2a In 
these calculations two approaches for the n — I image- 
state wave function have been used. First, it has been 
obtained as an eigenfunction of the model potential of 
Eq. (^), with the parameters chosen so as to reproduce 
the width and position of the energy gap and the binding 
energies at the f point (k|| — 0). Secondly, the n — 1 
image-state wave function has been obtained with use of 
the model potential of Eq. (^9|), but with the parameters 
chosen so as to reproduce the width and position of the 
energy gap and the binding energies at the correspond- 
ing values of k|| , thus allowing for the penetration of the 
image-state wave function into the bulk to increase with 
k|| . Though both approximations yield an image-state 
linewidth that increases with ky , it increases very slowly 
within the first approach and more rapidly within the 
second approach, showing the key role that the penetra- 
tion of the image-state wave function plays in the decay 
mechanism, i.e., as the coupling of the image-state wave 
function with bulk and intrinsic-surface states increases, 
the image-state linewidth is enhanced. 



As all theoretical calculations presented in Tables II, 
III and IV have been obtained with use of the wave func- 
tions of Eq. ( pO| ) , surface corrugation has not been taken 
into account. Estimating the influence of surface cor- 
rugation on the image-state broadening requires the use 
of three-dimensional wave functions for the evaluation of 
initial and final states and, also, for the evaluation of the 
screened interaction. Work along these lines is now in 
progress. 

An approximate way of including surface-corrugation 
effects on the final-state wave functions is to account for 
the actual effective mass of the intrinsic-surface and all 
bulk states, as obtained from the theoretically or experi- 
mentally determined dispersion of these states_In Table 
we compare the results of this calculationO^Ej with re- 
cent accurate TR-2PPE measurements of the lifetime of 
image states on Cu surfaces, showing very good agree- 
ment between theory and experiment. 

The present state of a theory that uses the self-energy 
formalism in combination with an accurate description, 
within a quasi one-dimensional model, of the key as- 
pects of image states, namely, the position and width 
of the energy gap and the binding energies of the intrin- 
sic and image-potential surface states has been shown to 
give quantitative account of the lifetime of image states 
on metal surfaces. Moreover, this theory also gives— a, 
linewidth of the Shockley surface state on Cu(lll)E£il 
at the F point that is in excellent agreement with re- 
cent very-highjXjSDintion angle-resolved photoemission 
measurements 



111 11; The calculated inelastic linewidth 



has been_£pund to be of 26 paeV, while measuremfints give 
30 mevH and 21±5 meV.Ela This calculationEi empha- 



sizes the extremely important role that the intrinsic sur- 
face state plays in the decaying mechanism of this state 
at the F point, resulting in a contribution of ~ 70% of the 
total linewidth. This surface-state contribution 
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aiams 
'^^^l and 



the difference between the experimental data 
theoretical results obtained within a bulk description of 
the broadening mechanism. 

If dielectric layers are grown on the metal substrate, 
one can analyze the layer growth by simply lookingj-at 
the energetics and lifetimes of image-state electrons.E^ 
Deposition of an overlayer on a metal substrate can 
change drastically the properties of image states, such 
as binding energies, wave functions, and lifetimes. This 
change depends on weather the adsorbate is a transition 
metal, an alkali metal, or a noble gas atom. In partic- 
ular, deposition of alkali metal adlayer on Cu(llJ-) de- 
creases the work function by nearly a factor of 2Ji^ The 
linewidths of image states on a single layer of Na and K 
on Cu(lll), Fe(llO) and Co(OOOl) were measured with 
2PPE by Fisher et aZjl£3 A large value of 150 meV was 
obtained for the first image state, which is quite close to 
the linewidth of the n = 1 image-state on Fe(llO) and 
Co(OOOl) but much larger than the linewidth of the n = 1 
image state on Cu(lll).U All these values were measured 
with 2PPE, and they include both energy and phase re- 
laxation contributions. Additionally, these experiments 
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showed the presence of the n = intrinsic surface state 
generated by the Na/K layer, which replaces the intrinsic 
surface state on Cu(lll). More accurate measurements 
of the intrinsic linewidth can be obtained with use of 
TR-2PPE spectroscopy. Nevertheless, the influence of 
impurities and imperfections on the linewidth remains to 
be evaluated. For more realistic estimates of the intrinsic 
linewidth, accurate models and/or first-principles calcu- 
lations arej_o£cessary. The same applies to other metal 
over layers. 

Overlayers of Xe and Kr on Ag(lll), Cu(lll), 



and 



163 



iutly with use 
spectroscopies. 



RumqOii have been 
of 2PPBiyil£a and TR2PPI 
All these measurements have shown that the lifetime of 
the n = 1 state increases significantly upon deposition of 
the noble atom adlayer on all metal substrates of inter- 
est. Qualitatively, this increase can be explained by the 
fact that the interaction of the image-state electron with 
the closed Xe or Kr valence-shell is repulsive and, there- 
fore, the probability amplitude of this state moves away 
from the crystal, as compared to the simple case of clean 
metal surfaces. Therefore, the coupling to the substrate 
decreases and the lifetime increases. The same qualita- 
tive argument can also explain the decrease of the binding 
energy of image states upon depositionpOj£l_X|e or Kr on 
metal substrates. Moreover, Harris et aZLSE£3 have stud- 
ied the evolution of image states as a function of the num- 
ber of deposited atomic layers of Xe on Ag(lll). They 
have found that with increase of a number of Xe layers 
the n — 2 and n = 3 image states evolve into quantum- 
well states of the overlayer. A qualitative interpretation 
of this behaviour of image states has been swen, within 
a macroscopic dielectric continuum model.ll£3 Unfortu- 
nately, no microscopic investigation of the image-state 
evolution of adlayers on metal surfaces that takes into 
account the band structure of both the substrate and 
the overlayer has been carried out. 

Defects on the surface or adsorbed particles cause elec- 
tron scattering processes that lead to phase relaxation 
of the wave function. This can be monitored in real 
time, in order to extract relevant information. In fact, 
measurements of the n — 1 and n = 2 image-potential 
states of CO adsorbed in Cu(lOO) indicate a decreasing 
dephasing time when the CO npnk'rii.lcs form an ordered 
c(2x2) structure on thc-S|Urface.EfZlE£3 Furthermore, mea- 
surements on Cu(lOO)E£0 have shown correlation between 
decay and dephasing, on the one hand, and the existence 
of surface defects, on the other hand. A first-principles 
description of this problem is still lacking, due to intrin- 
sic difficulties in dealing with the loss of two-dimensional 
translational symmetry. 

Another important field of research is the understand- 
ing of the processes leading to the electronic relaxation 
in magnetic matecials— The spin-split image states on 

offer the possibility of extract- 



magnetic surfacest£i 
ing information about the underlying surface magnetism. 
These spin-split states can decay in different ways and, 
therefore, their linewidths can be different. In particu- 



167 



lar, spiifetesolved inverse photoemission experiments on 
Fe(110)E£3 give an intrinsic linewidth of 140 (70) meV 
for the first minority (majority) image state. The differ- 
ence in the lifetime is of the order of the total linewidth 
of the n = 1 image state on other metal surfaces. At 
the same time, as the spin-splitting is only of ~ 8% of 
the total binding energy {Ei = — 0.73eV), it is unlikely 
that this splitting is responsible for the large difference 
between linewidths. Hence, one has to resort to details 
of the phase space of final states and to the screened 
Coulomb interaction as responsible for this effect. Work 
along these lines is now in progress. 

All these problems are of technological relevance and 
pose technical and theoretical questions that need to be 
answered in order to make a correct interpretation of 
what is really being measured. One technique is based on 
the ab initio description of the fast-dynamics of a^ave- 
packet of excited electrons in front of the surface.t£3 The 



time evolution will pick up all the relevant information 
concerning scattering processes and electronic excilations 
that can be mapped directly with experiments.Ej On a 
more complex and fundamental level, there is the theoret- 
ical description of coupled electron-ion dynamics, which 
is relevant in many experiments. 



VI. FUTURE 

We present here a brief summary of on-going and fu- 
ture work in the field of inelastic electron scattering in 
solids and, in particular, in the investigation of electron 
and hole inelastic lifetimes in bulk materials and low- 
dimensional structures. The advance in our knowledge 
is closely linked to the experimental developments that 
combine state-of-the-art angle-resolved 2PPE with ultra- 
fast laser technology. These investigations might be rele- 
vant for potential technological applications, such as the 
control of chemical reactions in surfaces and the devel- 
oping of new materials for opto-electronic devices. 

A theoretical and experimental challenge is the de- 
scription of the reactivity at surfaces. Experiments are 
being performed nowadays directed to get a deeper un- 
derstanding of the electronic processes involved. We note 
that electronic excitation is the initial step in a chemi- 
cal reaction, and the energetics and lifetimes of these 
processes directly govern the reaction probability. For 
example, we can achieve chemical selectivity thcoiigh a 
femtosecond activation of the chemical reaction.E£j This 



shows clearly that nonrandom dissociation exists in poly- 
atomic molecules on the femtosecond time-scale, by excit- 
ing the reactant to high energies (well above the thresh- 
old for dissociation) and sampling the products on time 
scales that are shorter than the rate for intramolecu- 
lar vibrational energy-distribution (this concept is rel- 
evant in chemical reactivity and assumes ergodicity or, 
equivalently, that the internal energy is statistically re- 
distributed). The idea of ergodicity has to be revisited 
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in this short time-scale. 

Very recently, it has been shown that selective adsorp- 
tion of low-energy electrons into an image-potential state, 
followed by inelastic scattering and desorption, can pro- 
vide information ou-the interaction between these states 
and the substrate.EZ3 A deep theoretical analysis of this 
interaction, as well as the role of the substrate/adsorbate 
band structure, is still lacking and is needed in order to 
interpret the experimental data. 

So far, we have concentrated our attention to the in- 
vestigation of bulk and surfaces, i.e., extended systems. 
From a technological point of view, and due to the fast 
miniaturization of the magneto- and opto-electronic com- 
ponents in current devices, the study of the electron dy- 
namics in nanostructures is of relevance. For example, 
alkali metals that have image states as resonances would 
have, in a finite piece of material, a well-defined image 
state with a long lifetime. These states are spatially lo- 
cated outside the nanostructure and, at least in princi- 
ple, could be used in a possible self-assembling mecha- 
nism to build controlled structures made of clusters, and 
also as an efficient external probe for chemical charac- 
terization. Measurements on negatively charged clusters 
would be able to assess this effect, as well as its size 
dependence. Experiments performed on a Na,^]^ clus- 
ter have looked at two decay mechanisms for the col- 
lective excitation, namely, electron and atom emission. 
Tha-estimated electronic escape time is of the order of 



1 fs.lilj The relaxation time for two-electron collisions in 
small sodium clusters has been estimated theoretically 
at the level of a timeijlependent local-density-functional 
approach (TDLDA)JlZj The computed values are in the 
range of 3 — 50 fs, which are between the direct elec- 
tron emission and the ionic motion (> 100 fs). These 
values compete with the scale for Landau damping (cou- 
pling of the collective excitation to neighbouring particle- 
hole states). A first non-perturbative approach to the 
quasiparticle l ifetim e in a quantum dot has been pre- 
sented in Ref. 173, where localized (quasiparticle states 
are single-particle-like states) and delocalized (superpo- 
sition of states) regimes are identified. Furthermore, 
if we wish to use these nanostructures in devices, we 
need to understand the scattering mechanism that con- 
trols the electronic transport at the nanoscale level. We 
expect new physical phenomena to appear in detailed 
time-resolved experiments in these systems, related to 
quantum confinement. In summary, the investigation of 
electron-electron interactions in nanostructures is still in 
its infancy, and much work is expected to be done in the 
near future. In particular, we are planning to investi- 
gate electron lifetimes in fuUerene-based materials, such 
as Ceo and carbon nanotubes. 

Asymmetries in electron lifetimes arise from the dif- 
ferent nature and localization of electrons; in this sense, 
noble and transition metals offer a valuable framework 
to deal with different type of electrons that present var- 
ious degrees of localization. New theoretical techniques 
should be able to address the calculation of excitations 



and inelastic electron lifetimes, including to some ex- 
tent electron-phonon couplings [which might be impor- 
tant and even dominant for high enough temperatures 
and very-low-energy electrons] and also both impurity 
and grain-boundary scattering. Final-state effects have 
been neglected in most practical implementations, and 
they might be important when there is strong localiza- 
tion, as in the case of transition and rare earth materials. 
In the case of semiconductors, electron-hole interactions 
(excitonic renormalization) strongly modify the single- 
particle optical absorption profile, asd they need to be 
included in the electronic response.lll3 Although similar 
interactions are expected to be present in metals, the 
large screening in these systems makes their contribution 
less striking as compare to the case of semiconductors. 
However, in the case of low-dimensional structures they 
might play an important role in the broadening mecha- 
nism of excited electrons and holes. 

All calculations presented in sections IV and V stop at 
the first iteration of the GW approximation. Although 
going beyond this approximation is possible, this has to 
be done with great care, since higher-order correctii; 
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tend to cancel out the effects of selfconsistencyt 
(see Appendix B). As we start from an RPA-like screen- 
ing, the net effect of including the so-called vertex correc- 
tions for screening electrons would be a reduction of the 
screening. Furthermore, a simpler and important effect 
to be included in the present calculations is related to 
the renormalization of the excitation spectral weight due 
to changes in the self-energy close to the Fermi surface. 
We know JJial—Lhis renormalization could be as_Large as 
0.5 for N£laO and of the order of 0.8 for Si.H This 
modifies the energy of the excitation and, therefore, the 
lifetime. We aim to include such effects in the calcula- 
tion of the inelastic eleetronic scattering process in noble 
and transition metalstSd, along the lines described in the 
Appendix B. The main idea is to work directly with the 
Green function in an imaginary-time/energy representa- 
tion. The choice of representing the time/energy depen- 
dence on the imaginary rather than the real axis allows 
us to deal with smooth, decaying quantities, which give 
faster convergence. Only after the full imaginary-energy 
dependence of the expectation values of the self-energy 
operator has been established do we use a fitted model 
function (whose sophistication may be increased as neces- 
sary with negligible expense) , which we then analytically 



continue to the real energy ax: 
citation spectra and lifetimes 



ider to compute ex- 
Furthermore, this 



technique is directly connected with finite-temperature 
many-body Green functions, and can be used to directly 
address temperature effects on the lifetime that can be 
measured experimentally. 

An interesting aspect in the theory of inelastic elec- 
tron scattering appears when one looks at the energy 
dependence of the electron lifetime in layered materi- 
als. In a semimetal as graphite, the lifetime has been 
found to be irpteisely proportional to the energy above 
the Fermi levelt®^, in contrast to the quadratic behaviour 
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predicted for metals with the use of Fermi-hquid theory 
(see Section IV). This behaviour has been interpreted in 
terms of|-electron-plasmon interaction in a layered elec- 
tron gasjl23 however, this is not consistent with the fact 
that a laiceied Fermi-liquid shows conventional electron 



lifetimes.E£D A different interpretation based on the par- 



ticular band structure of graphite (with a nearly point- 
like Fermi surface) yielding a reduction of the screetu. 
ing can explain the linear dependence of the lifetime. a 
A similar linear dependence of the inelastic lifetime has 
been foimd, for other semiconducting-layered compounds 



as SnS2.li^ We are presently working on the evaluation, 



within the GW approxiniatjon, of electron lifetimes in 



these layered compounds .tiHI The special band structure 



of graphite has also been invoked as responsible for the 
peculiar plasmon dispersion and damping of the surface 
plasmon.E£il Therefore, a careful analysis of the layer- 
layer interaction and broadening of the Fermi surface 
needs to be included, in order to understand this be- 
haviour. We note that in a metal like Ni the imaginary 
part of the self-energy shows a quadratic Femii-|liquid be- 
haviour, which becomes linear very quicklyt 
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Together with the self-energy approach discussed in 
this review, an alternative way of computing the excita- 
tion spectra of a many-body system, which is based on 
information gleaned from an ordinary ground-state cal- 
culation, ifjJJiR-tiine-dependent density-functional theory 



(TDDFT)Jl£3-Ea In this approach, one studies how the 
system behaves under an external perturbation. The re- 
sponse of the system is directly related to the A^-particle 
excited states of an iV-particle system, in a similar man- 
ner as the one-particle Green function is related to the 
{N + 1)- and (TV — l)-particle excited states of the same 
system. TDDFT is an ideal tool for studying the dy- 
namics of many-particle systems, which is based on a 
complete representation of the XC kernel, K^'^, in time 
and space^ — .One computes the time-evolution of the 
systemE23T23 without resorting to perturbation theory 
and dealing, therefore, with an external field of arbitrary 
strength. The fact that the evolution of the wave func- 
tion is mapped for a given time-interval helps one to ex- 
tract useful information on the dynamics and electron 
relaxation of many-electron systems. The method does 
not stop on the linear response and includes, in princi- 
ple, higher-order nonlinear response as well as multiple 
absorption and emission processes. 

On a more pure theoretical framework, the connection 
between TDDFT and many-body perturbation theory is 
needed, in order to get further insight into the form of 
the frequency-dependent and non-local XC kernel K^'^. 
If one were able to design an XC kernel that works for 
excitations as the LDA does for ground-state properties, 
then one could handle many interesting problems that are 
related to electron dynamics of many-electron systems. 

In summary, many theoretical and experimental chal- 
lenges related to the investigation of lifetimes of low- 
energy electrons in metals and semiconductors are open, 
and even more striking theoretical and experimental ad- 



vances are ready to come in the near future. Lifetime 
measurements can be complementary to current spectro- 
scopies for the attainment of information about general 
properties (structural, electronic, dynamical, ....) of a 
given system. 



ACKNOWLEDGMENTS 

The authors would like to thank I. Campillo, M. A. 
Cazalilla, J. Osma, I. Sarria, V. M. Silkin, and E. Zarate, 
for their contributions to some of the results that are re- 
ported here, and M. Aeschlimann, Th. Fauster, U. Hofer, 
and M. Wolf, for enjoyable discussions. Partial support 
by the University of the Basque Country, the Basque 
Unibertsitate eta Ikerketa Saila, the Spanish Ministerio 
de Educacion y Cultura, and Iberdrola S. A. is gratefully 
acknowledged. 



APPENDIX A: LINEAR RESPONSE 

Take a system of interacting electrons exposed to 
an external potential V'^^'^{r,uj). According to time- 
dependent perturbation theory and keeping only terms 
of first order in the external perturbation V'^^*{r,uj), the 
charge density induced in the electronic system is found 
to be 



o"^(r,tj) 



dr'x(r, r';w)F'=^*(r',t^), 



(Al) 



where x(r, i"'; 'j-') represents the so-called linear density 
response function 



X(r,r',^) = ft*o(i")Pno(r') 

n 

1 



1 



LO - LOnO + W 



(A2) 



Here, Uno — En — Eq and p„o(r) represent matrix el- 
ements taken between the unperturbed many-particle 
ground state |^'o) of energy Eq and the unperturbed 
many-particle excited state of energy En'. 



p„o(r) = (*„|p(r)|*o), 
p(r) being the charge-density operator, 

N 



p(r) 



(A3) 



(A4) 



i=l 



and Ti describing electron coordinates. 

In a time-dependent Hartree or random-phase approx- 
imation, the electron density induced by the external po- 
tential, F'^^*(r,a;), is approximated by the electron den- 
sity induced in a noninteracting electron gas by the total 
field y^^*(r,w) + V"™'^(r,w): 



16 



V""^{v',uj)\. (A5) 

This approximation for the induced electron density can 
be written in the form of Eq. ( Al), with 



X 



RPA 



(r,r';w) = x"(i":i"';^) + ^dri j Ar2 

xx°(r, ri; a;) v{v, ~ v,) x^^-'lrs, r'; cu), (A6) 



where x'^(r, r'; lu) is the density-response fmiction of non- 
interacting electrons, 



X°(r,r';c.) = 2^ 



Si - £j + (uJ + iri) 



x^,(r)^*(r)^,(r')^*(r'), (A7) 
(pi (r) representing a set of single-particle states of energy 

In th&iramework of time-dependent density-functional 
theory,li2j~t2j the theorems of which.,aen,eralize those of 



the usual density-functional theory. 
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the density re- 



sponse function satisfies the integral equation 

x{r,r';u;)^X°ir,r';cj) + J dri J dra x°(r, ri; w) 

X [«(ri-r2)+if^^(ri,r2;w)]x(r2,r';a;), (A8) 

the kernel K^''{r, r'; lu) representing the reduction in the 
e-e interaction due to the existence of short-range^JLC 
effects. In the static hmit (lu ^ 0), DFT shows that££3 



X^^(r, r';Lu^ 0) 



6n{r)6n{r') 



(A9) 



no(r) 



where £'a;c['T'] represents the XC energy functional and 
no(r) is the actual density of the electron system. 

In the case of a homogeneous electron gas, one intro- 
duces Fourier transforms and writes 



P. 



q,a; 



q,uj • 



Within RPA, 



RPA ^ , RPA 
Aq.uj Aq.cj ' Aq,a; QAq.cj i 



(AlO) 
(All) 



where 



1 
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^k+q - Ek + IT] 



(A12) 



Vq — Air/q'^ is the Fourier transform of the Coulomb po- 
tential, Ek = and rik are occupation numbers, as 
given by Eq. (|). 

In the more general scenario of TDDFT, 



Xq,' 



Aq,cj 
Aq.o; 



Xq,' 



Xq,c 



Gci.oj) Xq t 



(A13) 



where K^"^^ is the Fourier transform of the XC kernel 



K'"'{r,r',uj). InEq. ([A1^), we have set 



^q Gq,a;i 



(A14) 



Gq,(j being the so-called local-field factor. In the local- 
density approximation, which is rigorous in the long- 
wavelength limit {q ^ 0), it follows from Eqs. (A9) and 
(|A14|) that 



where 



f^LDA _ 4 I -J 



^ _ 1 Att (fE^ 
4 q% dnl 



(A15) 



(A16) 



Ecina) being the correlation contribution to the ground- 
state energy of the uniform electron gas. This quan- 
tity has been extensively studied, rangip g from the 
simple Wigner intemDlation formulaoE23 to accurate 
parametrization£2il'E23 based— on Monte Carlo calcula- 



tions by Ceperley and Alder 

Diffusion Monte Carlo calculations of the static 
dei^aLtv^ixesponse function of the uniform electron 



gast 



y-T'e« 

: 204|g05| i^ave shown that the LDA static local-field fac- 



tor of Eq. ( Al5| ) correctly reproduces the static response 
for all (J < 2 gi?, as long as the exact correlation energy 
is used to calculate A. For larger values of q both exact 
and RPA density- response functions decay as l/(j^, their 
difference being of order l/g**, and fpe, details of Gq.o are 
expected to be of little importance£231 
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Calculations of the frequency dependence of the looak 
field factor were carried out by Brosens and coworkers^ 
and, more recently, by Richardson and Ashcroft.tiH 
Local-field factors are in general expected to be complex 
for nonzero frequencies, and the importance of their fre- 
quency dependence is refiected, e.g., in the finite lifetime 
of the volume plasmon. For aluminum, measurements of 
Gq,w have shownHJ a weak uj dependence of the local- 
field factor for energies below ~ 35 eV. 

If the homogeneous electron gas is exposed to an ex- 
ternal test charge of density /9^^'(r, t), one writes 



q,a; 



q ^q,cj' 



(A17) 



and with the aid of Eq. ( |AlCl| ) one finds the Fourier trans- 
form of the total change of the charge density /3*°*(r, t) = 
pext^j.^ i) + p™'^(r, t) to be given by the following expres- 
sion: 



-1 ext 
q,cj r'q,:^ ' 



(A18) 



where eq ^^ is the so-called test_charge-test_charge dielec- 
tric function: 
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q,CL; 



1 



(A19) 



Hence, this dielectric function screens the potential both 
generated and 'felt' by a test charge. 

If the external potential is that generated by an elec- 
tron, then one writes 



K?rJ=«q(l-Gq,.)p: 



(A20) 



and with the aid of Eq. (|A10D one finds Eq. (|A1S| ), 
but iiow—£q,uj being the test_charge-electron dielectric 
functionO 



1 + Wq(l - Gq,u;)XqA 



(A21) 



By combining Eq. (A13) with either Eq. ( Aig| ) or Eq. 
(A21), one can easily obtain the following expressions for 
the test_charge-test_charge and the test_charge-electron 
dielectric functions. 



and 



eq,L = 1 



1 + Wq Gq,^ Xq,i 



Wq(l - Gq,^)Xq,, 



(A22) 



(A23) 



respectively. If t he lo cal-fi eld fa ctor G^^^j is set equal to 
zero, bothjEfls. ( [A23| ) and ( [A24| ) yield the RPA dielectric 
functionESEa 



RPA 
q,aj 



1 



(A24) 



In terms of e^^^, Eqs. ( |A23|) and (|A24|) yield Eqs. (|E 
and (p^), respectively. 



APPENDIX B: GW APPROXIMATION AND 
BEYOND 



Let us introduce the many-body Green function:El^ 



Gir t, r' t') - -leit - t') «|^(r, t)^^ (r', 

+\6{t'-t) (vl/^|V/(r',i')V^(r,i)lO. (Bl) 

In this equation, ■0'l'(r, stands for a (A^ -I- 1)- 

electron state in which an electron has been added to the 
system at point r and time t. When t' < t, the many- 
body Green function gives the probability amplitude to 
detect an electron at point r and time t when a (possibly 
different) electron has been added to the system at point 
r' and time t'. When t' > t, the Green hmction describes 
the propagation of a many-body state in which one elec- 
tron has been removed at point r and time t, that is, the 
propagation of a hole. 

For a system of interacting electrons, there is little 
hope of calculating Q{r,r',uj) exactly. One usually has 
to resort to perturbation theory, starting from a suitably 
chosen one-electron problem with a Hamiltonian iJo(r), 



eigenfunctions 0i(r), and eigenenergies Si. Hence, the 
noninteracting Green-iunction t/'^(r, r', u) is given by the 
following expression :t22l 



uj - et + hj sgn(ei - Ep) 



(B2) 



In usual practice the LDaEiI is considered, which pro- 
vides a local one-electron potential, ulda{t^)- 

The exjadj Green function obeys the following Dyson's 
equationJ123 



(B3) 



^ + ^V2-u(r)jg(r,r',c^) + 
j dr"S(r, r",Lu)g{r", r, lu) = 5(r - r'). 



where the integral kernel I](r, r', w) is known as the self- 
energy. It can be understood as the complex non-local 
energy-dependent potential felt by the electron added to 
the system at point r' and time t'. This potential arises 
from the response of the rest of electrons to the presence 
of the additional electron. However, one must be care- 
ful with this interpretation, since the many-body Green 
function for t' < t not only describes the propagation 
of the additional electron, but also that of the whole 
(A^ -I- l)-electron system. This means, for instance, that 
the self-energy also accounts for the exchange processes 
that occur in a system of indistinguishable particles. 

To obtain the inelastic lifetime of one-clectron-like ex- 
citations, called quasiparticles, one seeks for the poles of 
the many body Green function. A good estimate of the 
position of these poles can be obtained by projecting Eq. 
(B3) onto the chosen basis of one-electron orbitals, and 
neglecting the off-diagonal terms in the self-energy, i.e.. 



e,-AI],,H«0, 



(B4) 



where 



AE,,H = JdrJ dr'0*(r) [S(r,r';a;) 

~ULDA{r)S{r~r')]Mr')- (B5) 



In general, Eq.(B4) has complex solutions ai uj ~ 
E — ir/2, where E is the excitation energy and F the 
linewidth of the quasiparticle. Near the energy-shell 
{uj ^ £i, Ei being the free-particle energy), one finds for 
the linewidth r^^ ~ F: 



where 



Z,, = 



1 - 



aRcAS,,(w) 



(B6) 



(B7) 



is a renormalization constant. On the energy-shell {Zi 
1), one writes 
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-2ImS,,(e,), 



(B8) 



and after noting the reahty of the matrix elements of the 
LDA potential one finds Eq. (p^. j— . 

Within many-body perturbation theory,E3 it is possible 
to obtain S(r,r',a;) as a series in the Coulomb interac- 
tion v{r — r'). Due to the long range of this interaction, 
such a perturbation series is expected to contain diver- 
gent terms. However, it has been known for a long time 
that when the polarization induced in the system by the 
added electron is taken into account the series is free 
of divergences. Thus, the perturbation series for the self- 
energy can be rewritten in terms of the so-called screened 
interaction W{r,r';Lu). 

To lowest order in the screened interaction, the self- 
energy reads: 

S(r,r',c.) = i J ^e'^'^'g{r,r',u; + u;')Wir,r',u;), 

(B9) 

which is the so-called GW approximation. The screened 
interaction is given by Eq. (p^), in terms of the ex- 
act time-ordered density-response function of interacting 
electrons, or, equivalently, by the following integral equa- 
tion 

W{r, r',uj) = v{r — r') + 
J dri J dr2 «(r-ri)n(ri,r2,cj)VK(r2,r',u;), (BIO) 

where n(r, r',a;) represents the polarization propagator. 
In the GW approximation, 

nG^(r,r',c.) = -i J ^giry,Lu')gir\r,u;'-Lu). 

(Bll) 

If Q is replaced in this expression by Q'^ (n*^*^) ^ II"), 
one easily finds 



and 



Ren"(r,r';w) =Rex°(r,r';w) (B12) 



Imn°(r,r';w) = sgnwlmx"(r,r'; w), (B13) 



where x°(i": i"'; represents the retarded density- 
response func tion of noninteracting electrons, as defined 
by Eq. (A7). For positive frequencies {uj > 0), both 
n°(r,r';tj) and x°(r, r';w) coincide. 

The GW approximation gives a comparatively simple 
expression for the self-energy operator, which allows the 
Green function of an interacting many-electron system 
to be computed by simply starting from the Green func- 
tion Q*^ {r , r' ; uj) of a fictitious system with an effective 
one-electron potential. The GW approximation has been 
shown to be physically well motivated, specially for met- 
als where the Hartree-Fock approximation leads to un- 
physical results. 



Eqs. (B3)-(B11) form a set of equations which must 
be solved self-consistently for Q{r,r',uj). This means 
that the Green function used to calculate the self-energy 
must be found to coincide with the Green function ob- 
tained from the Dyson equation with the vet^L-same self- 
energy. However, there is some evidenca]£d support- 
ing the idea that introducing the noni nteracting Green 
function ^°(r,r';w) both in Eq. (^) and Eq. ( |bTi1 ) 
{G'^W^ approximation) one obtains accurate results for 
the description of one-electron properties such as the ex- 
citation energy and the quasiparticle lifetime. However, 
self-consistency modifies the one-electroii-escitation spec- 
trum (excitation energies and lifetimes) Jilj'EZZI as well as 
the calculated screening properties. Self-consistent calcu- 
lations have been performed only very recently for the ho- 
mogeneoois electron gaslll3, simple semiconductors, and 
metals 
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Discrepancies between G'~'W'~' and self-consistent GW 
calculations seems to be origin.atftd in the fact that the 
so-called vertex corrections Ji2^Hi2d which go beyond the 
GW approximation, need to be included as well. In- 
clusion of these corrections might cancel out the effect 
of self-consistency, thereby full self-consistent self-energy 
calculations yielding results that would be close to G°W° 
results. The main outcome of self-consistent GW-calcula- 
tions for the electron gas is that the total energyl£i3 turns 
out to be strikingly close to the total energy calculated 
with use of quantum Monte Carlo techniques E23 This 
result may be related to the fact that the self-consistent 
GW scheme conserves electron-number, energy, and total 
momentum, that is, fulfills the microscopic conservation 
laws.Eil 

The simplest improvement to the GW approximation 
is achieved by introducing a vertex correction that is 
consistcBt-with an LDA calculation of the one-electron 
orbitals,ci3 the XC potential being regarded as a self- 
energy correction to the Haxtree approximation. Based 



on this idea, the vertex YrlB can be easily expressed in 



terms of the static local field correction of 
This is the so-called GWF approximation£23TLil3 In 
this approximation, the polarization propagator is for- 
mally equivalent to the density-response function ob- 
tained within linear response theory in th^Jt^piejWork of 

In the 
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time-dependent density-functional theory, 
case of a homogeneous electron gas, this approximation 
yields the test_charge-electron dielctric function of Eq. 



APPENDIX C: LIFETIMES OF HOT 
ELECTRONS NEAR THE FERMI LEVEL: 
APPROXIMATIONS 

The damping rate of hot electrons near the Fermi level 
{E « E~ Ef) is obtained, within RPA, from Eq. ^ 
with the inverse dielectric function of Eq. (|35|), i.e.: 
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dq , ,_2 {E - EfY 
7''^'"' h ■ 



For small values q {q « 2qp) 



with 



(CI) 



(C2) 



(C3) 



q/2qp, and (3 — 1/3, and one obtains 



/— -3/2 



8^1 -/3c 



where 



ao 



1 

arctan — 
ao 



a 



1 + 



{E - Ef) 



(C4) 



(C5) 



This is the expressipH first obtained by RitchieEj and by 
Ritchie and Ashley.EJ 

In the high-densi ty li mit {qp — > 00), the static dielec- 
tric function of Eq. ( C2 ) is 



TF 
^q,0 



1 



Itf 



(C6) 



which can also be derived within the Thomas-Fermi sta- 
tistical model, qTF — \f^^FT^ be ing the Thom as-F ermi 
momentum. By introducing Eq. ( |C6| ) intOjEp. (CI) one 
obtains the expression derived by Quinn,lHJ whi ch can 
also be obtained by just taking /3 = in Eq. (|C4|): 



r- -3/2 



arctan — 
a 



{E-Epf 



(C7) 



If one replaces the static dielectric function of Eq. (|C2| ) 
by t he high-density limit [qp —^ 00), as given by Eq. 
(C6), and extends, at the same time, the integration of 
Eq. (CI) to infinity, one finds Eq. (36), which can also 
be obtained by just keeping the first-order term in the 
expansion of Eq. {C4) in q^^ . If we further replace ki — > 
qp in Eq_J3q), then the formula of Quinn and Ferrell is 
obtained J2a 



with 



TQ'p = Cir,)iE~Ep)' 



or, equivalently, 



128 Ej,'' 



Cir ) (llH^r^n 



(C8) 



(C9) 



(CIO) 



In Fig. 18, the ratio t/tqp is plotted against for hot 
electrons in the immediate vicinity of the Fermi surface 
{ki ~ qp), as computed from either Eq. ( pT] ) or Eq. (CI) 
(solid line), and also from Eqs. (C4) (dotted line) and 
( |C7| ) (dashed-dotted hne). In the range < < 6, 
Eq. (C4) reproduces the full RPA calculation within a 
2%, while either Eq. (36) or Eq. ( |C^ ) reproduce the 
full calculation within a_.7%. Differences between the 
approximation of QuinnEj (Eq. (|C7|)) and the full RPA 
calculation go up to 14% at r^. = 6. 

Departures of the predictions of Eq. (CS) (dotted 
line) from the full RPA calculation (solid lineTarise from 
small differences between the dielectric function of Eq. 
and the actual static RPA dielectric function. Fur- 



ther replacing the dielectric function of Eq. ( p2| ) by the 
high-density limit (qp — > 00), as given by Eq. ( |C(^ ), 
leads to a too-strong Thomas-Fermi-like screening and 
results, therefore, in a large overestimation of the life- 
time (dashed-dotted line). However, this overestimation 
is largely compensat ed if one also takes gj;^ ^ 00 in the 
integration of Eq. (CI), thereby the lifetime of Quinn 
and Ferrell tqp (dashed line) being closer to the full 
RPA calculation than the lifetime derived from Eq. (C7) 
(dashed-dotted line). 
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TABLE I. Available experimental data for hot-electron 
lifetimes in metals, as obtained by time-resolved two-photon 
photoemission and ballistic electron emission microscopy 
(BEEM). 



Metal 


Reference (technique) 


Cu 


[ 17|,[ 1^, 2C1 (TR-TPPE) 


Ag 


[ 17|,[ |l6[ (TR-TPPE) 


Au 


[ 17j,[ 2| (TR-TPPE); [ 23\ (BEEM) 


Ta 


[ 17|,[ L6| (TR-TPPE) 


Pd 


[Q4\ (BEEM) 


Al 


[ 25 1 (TR-TPPE) 


Co 


[ 2eW 27] (TR-TPPE) 


Fe 


[0 (TR-TPPE) 



TABLE IL Linewidth (inverse lifetime) of image states, in 
meV. 



Surface 


Image state 


2PPE 


TR2PPE 


Theory 


Cu(lOO) 


n=l 


28 ±6" 


16.5 ± 3/2"''= 
5.5 ± 0.8/0.6''''= 


18"; 26''; 22"= 




n=2 




5"= 




n=3 




2.20±0.16/0.14''''= 


1.8"= 


Cu(lll) 


n=l 


16±4-';85± lO^'S 


38 ± 14/9''; 30' 


20-';421"; 118''; 38"= 


Ag(lOO) 


n=l 


21 ±4" 


26 ±18/7*=; 12 ±1'= 


22"; 25'' 




n=2 


3.7 ±0.4" 


3.7 ±0.4'-'; 4.1 ±0.3/0.2'=; 


5" 




n=3 




1.83 ±0.08'= 




Ag(lll) 


n=l 


45 ±10"; 55' 


22 ± lo/e*" 


58'; 123"; 110'' 


Au(lll) 


n=l 


160 ± 40" 




617" 


Pd(lll) 


n=l 


70 ± 20" 




40"; 35" 


Ni(lOO) 


n=l 


70 ±8" 




24" 


Ni(lll) 


n=l 


84 ± 10" 




40" 


Co(OOOl) 


n=l 


95 ± 10" 




40" 


Fe(llO) 


n=l 


130 ± 30" 




95" 
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^Ref. 33 
=Ref. 35 
Bl 
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•^Ref. 
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TABLE III. Calculated contributions to the linewidth, in 
meV, of the n — 1 image state on Cu surfaces. 



Surface 






Winter 


r 


Cu(lOO) 


24 


14 


-16 


22 


Cu(lll) 


44 


47 


-54 


37 



TABLE IV. Linewidths Ti and La, in meV, of the n = 1 
image state on Cu(lll), as calculated for non-zero momenta 
parallel to the surface and with use of two models for the 
n — 1 image-state wave function (see text and Ref. |7q). 



k (a.u.) 


Li 




0.0570 


38.1 


38.9 


0.0912 


38.5 


43.6 


0.1026 


38.7 


47.0 



TABLE V. Lifetimes of image states on Cu surfaces, in fs. 



Surface 


Image state 


TR2PPE 


Model potential 
calculation 


Cu(lOO) 


n=l 


40 ± 6"'" 


30^=; 38" 




n=2 


120 ± 15"''' 


132^^; 168'' 




n=3 


300 ±20"''' 


367'=;480'' 


Cu(lll) 


n=l 


18 ±5"=; 22 ±5-^ 


17.5'=;22.5'''S 



. 33, U. Hofer et al. 

35, L L. Shumay et al. 

77, E. V. Chulkov et al. 
. 30, L Sarria et al. 

30, M. Wolf et al. 
. 31, E. Knoesel et al, at low temperature, T 
. 78, J. Osma et al. 
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FIG. 1. Scattering of an excited electron with the Fermi 
sea. The probe electron is scattered from a state (j>i{r) of 
energy Ei to some other state 4>f{jc) of energy Ef, by carrying 
one electron of the Fermi sea from an initial state (j>i' (f) of 
energy E^i to a final state <j>f'{r) of energy Efi, according 
to a dynamic screened interaction W{r — r' , Ei — Ef). Ef 
represents the Fermi level. 



FIG. 7. Averaged lifetimes of hot electrons in Cu, versus 
E — Ef, as obtained from Eq. (^9|) with the local electron den- 
sity of Ref. ^ (dotted line), and from Eq. ( [ill ) with the model 
dielectricr^Linction of Eq. (kol) and the recipe described by Sal- 
vat et a/E3 to obtain the optical energy-loss function (dashed 
line). The solid line represents the result obtained from Eq. 
( Pt[ ) with use of the free-electron gas RPA energy-loss function 
and Ts = 2.67. 



FIG. 2. Ratio of the lifetime of electrons above the Fermi 
level {E > Ef) to the lifetime of holes below the Fermi level 
{E < Ef), as a function of |i5 — Ef\, calculated within RPA 
for an electron density equal to that of valence (4s^) electrons 
in copper, i.e., = 2.67. 



FIG. 3. Ratio t/tqf between the lifetime r evaluated in 
various approximations and the lifetime tqf of Eq. (^7|), 
versus E — Ef, as obtained for hot electrons in a homoge- 
neous electron gas with Vs — 2.67. The solid line represents 
the result obtained from Eq. (|ll|), within RPA. Results ob- 
tained from Eqs. (C4) and ( jp'/\ ) are represented by dotted and 

The dashed line represents 



dashed-dotted lines, respective^ 
the result obtained from Eq. (B6|) 



FIG. 4. Exchange and correlation effects on the lifetime of 
hot electrons with E — Ef = 1 eV. The dashed line repre- 
sents, as a function of rs, the ratio between lifetimes derived 
from Eq. ( |ll| ) with use of the dielectric function of Eq. ( [l^ ) 
with (Gq,(j ^ 0) and without (Gq,cj — 0) local-field correc- 
tions. The dotted line represents, as a function of Vs, the 
ratio between lifetimes derived from Eq. (|l^) with use of the 
dielectric function of Eq. (|l^) with {Gq,Lj 7^ 0) and with- 
out (Gq.ui = 0) local-field corrections. If the local-field factor 
Gq,u is taken to be zero, both Eqs. (|l|) and (|l|) give the 
same result (solid line). 



FIG. 5. As in Fig. 4, but for hot electrons in a homoge- 
neous electron gas with Vs = 2.67, and as a function of the 
electron energy E — Ef with respect to the Fermi level. 



FIG. 6. Averaged lifetimes of hot electrons in Al, versus 
E — Ef, as obtained from Eq. (|9|) with the local electron den- 
sity of Ref. ^ (dotted line) , and from Eq. ( [lH ) with the model 
dielectrio-imction of Eq. mm and the recipe described by Sal- 
vat et aliB to obtain the optical energy-loss function (dashed 
line). The solid line represents the result obtained from Eq. 
( pr] ) with use of the free-electron gas RPA energy-loss function 
and rs = 2.07. 



FIG. 8. Averaged lifetime-widths of excited electrons 
in Cu, versus E — Ef, as obtained from Eq. ( psj ) with the local 
electron density of Ref. jsi] (dotted line), and from Eq. ( |ll| ) 
with the model dielectric rfunction of Eq. (^o[) and the recipe 
described by Salvat et a/Ej to obtain the optical energy-loss 
function (dashed line). The dashed-dotted-dotted-dotted line 
represents the result obtained from Eg. ( prj) with use of 
the model dielectric function of Eq. ([4C|) and the experi- 
mental optical energy-loss function of Ref. |llO| . The dotted 
line represents the result of using the 'universal' reJejtionship 
r"^ = 0.13 {E - Ef) proposed by Goldmann et alH 



FIG. 9. Hot-electron lifetimes in Al. Solid circles represent 
the full ab initio calculation of t{E), as obtained after av- 
eraging T of Eq. (j^j over wave vectors and over the band 
structure for each wave vector. The solid line represents the 
lifetime of hot electrons in a FEG with Vs — 2.07, as obtained 
from Eq. (|ll|). The dotted line represents the statistically 
averaged lifetime, as obtpiped from Eq. (js^) by following the 
procedure of Tung et al !t3 



FIG. 10. Hot-electron lifetimes in Cu. Solid circles rep- 
resent the full ab initio calculation of r{E), as obtained after 
averaging r of Eq. (^^ over wave vectors and over the band 
structure for each wave vector. The solid line represents the 
lifetime of hot electrons in a FEG with rs — 2.67, as obtained 
from Eq. (|ll|). The dotted line represents the statistically 
averaged lifetime, as obtamed from Eq. (|3^) by following the 
procedure of Tung et a/.uOpen circles represent the result 
obtained from Eq. (|^ by replacing hot-electron initial and 
final states in |_Bi/(q + G)|^ by plane waves and the dielectric 
function in |eG,G(q, <i')|~^ by that of a FEG with = 2.67, 
but with full inclusion of the band structure in the calculation 
of Im [eG.G(q, ^^)]- Full triangles represent the result obtained 
from Eq. (p3) by replacing hot-electron initial and final states 
in \Bif{q + G)f by plane waves, but with full inclusion of the 
band structure in the evaluation of both Im [eG,G(q, t*^)] and 
|eG,G(q, 1.^)1"^. 
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FIG. 11. Scaled hot-electron lifetimes in Cu. Solid circles 
represent the full ab initio calculation of r{E), as obtained 
after averaging r of Eq. ( |2^ ) over wave vectors and over the 
band structure for each wave vector. The solid line represents 
the lifetime of hot electrons in a FEG with Vs — 2.67, as ob- 
tained from Eq. (pi]). The dashed-dotted line represents the 
statistically averaged lifetime, as obtakied from Eq. ( ^ ) by 
following the procedure of Tung et al £j The dashed line rep- 
resents the result of following the procedure described in Ref. 
and the dotted line is the result of using the 'universal' 
onship = 0.13 {E — Ep) proposed by Goldmann et 
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FIG. 18. Ratio t/tqf between the lifetime r evaluated in 
various approximations and the lifetime tqf of Eq. (^7[), 
versus the electron-density parameter r^, as obtained for hot 
electrons in the vicinity of the Fermi surface {E ^ Ef)- The 
solid line represents the result obtained, within RPA, from 
either Eq. Q or Eq. (^. Results obtained from Eqs. (|c|) 
and ( |C7[ ) are represented by dotted and dashed-dotted lines, 
respectively. The dashed line represents the result obtained 
from Eq. M). 



FIG. 12. Experimental lifetiisps of low-energy electrons in 
Cu, pj taken from Knoesel et alci (solid circles), from Ogawa 
et alc3 (Cu[100]: open circles, Cu[110]:i-apen squares, Cu[lll]: 
solid squares), from Aeschlimann et alUl (solid triangles), and 
from Cao et ali3 with oj/ot. = 1.63eV (open diamonds) 



FIG. 13. Schematic plot of the model potential of Eq. 
191). Vertical solid lines represent the position of atomic lay- 



FIG. 14. The probability amplitude of the n — 1 image 
state on (a) the (110) surface of Li, as obtained from the model 
potential of Eq. (|4S|) (solid line) and from pseudopotential 
calculations (dotted line), and (b) the (100) surface of Cu, as 
obtained from the model potential of Eq. (Us) (solid line) and 
from linear augmented plane- wave calculations (dotted line). 
Vertical solid lines represent the position of atomic layers. 



FIG. 15. Imaginary part of the electron self-energy, versus 
z, for three fixed values of z' (solid circles), as calculated for 
the n — 1 image state on (a) Cu(lll) and (b) Cu(lOO). 



FIG. 16. (a) Maximum value of Im [-E(z, z'; ky = 0; Ei)] 
(see the text) for the n = 1 image state on Cu(lll) (solid 
line) and Cu(lOO) (dotted line). Vertical lines represent 
the position of atomic layers in Cu(lll) (solid lines) and 
Cu(lOO) (dotted lines). (b) As in (a), for the separate 
contributions to the n = 1 image state on Cu(lll). The 
solid line with circles describes the total maximum value of 
Im S(z, z'; k|| = 0; -Ei)] , the dashed line represents the con- 
tribution coming from the decay into the intrinsic surface 
state, and the dotted line, the contribution from the decay 
into bulk states. The solid line represents the result of replac- 
ing the realistic model-potential 4>f{z) final wave functions en- 
tering Eq. ( p^ by the self-consistent jellium LDA eigenfunc- 
tions of the one-electron Kohn-Sham hamiltonian but with 
the restriction that only final states with energy Sf lying be- 
low the projected band gap are allowed. 



FIG. 17. Schematical representation of the electronic 
structure of Cu(lll) and Cu(lOO). 



28 



rigure i 




1 1 1 1 1 1 1 1 1 1 

: E-Ep=1eV 


1 1 1 1 1 1 1 1 

Eq. (13) 


1 1 1 1 


1 1 














RPA 






^ ^ ^q-(12) 






~ 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 









Figure 5 

1 ,2 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I r 



1,1 



r =2.67 

s 



Eq. (13) 



RPA 



0,9 



0,8 



- - - - Eg. (12) 



0,7 



0,5 1 1,5 2 2,5 3 3,5 

E - Ep (eV) 



Figure 7 



Cu 




1 1,5 2 2,5 3 3,5 4 

E- (eV) 



Figure 10 




0,5 1 1,5 2 2,5 3 3,5 4 

E - (eV) 



100 



> 

CD 

3 60 



C\J 

■'Hi 

LU 40 

e 20 




Figure 1 1 



1 1 1 1 1 1 

Cu 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

Eq.(39)__ - - 




• 


... Eq.(28) 




Zarate 








__FEG^ — -^^^^ 




1 1 1 1 


Goldmann 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 


1 1,5 


2 2,5 3 3,5 

E- (eV) 





> 



160 
140 
120 
100 
80 
60 
40 
20 




Figure 12 



□ 

□ □ 



□ 



o 



o 



o 



o 



8 



o 



o 



oo 



O 
O O 



□ 



o 



o 



o 



o 



•ooo 



^ □ □ □□ 



▲ ▲ 



Cu 



0,5 1 



1,5 



2,5 



3,5 4 



E-E^ (eV) 



Figure 14 



0.06 




z (a.u.) 



Figure 15 




Figure 16 

350 




-30 -25 -20 -15 -10 -5 5 10 



Z (a.u.) 



Figure 17 




